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ABSTRACT 

The cascade algorithm that is used for extended surface analysis depends on a 
new parameterization called the thermal transmission matrix to represent a single fin. 
This thermal transmission matrix, which is intended to replace the more familiar fin 
efficiency as a design and analysis parameterization, is a linear transformation that 
maps conditions of heat flow and temperature at the fin tip to heat flow and 
temperature conditions at the fin base. The cascade algorithm was derived by 
resorting to an analogy between a fin and the electrical transmission line. The 
cascade algorithm permits a fin to be subdivided into many subfins each having a 
thermal transmission matrix and then the individual transmission matrices for each of 
the subfins can be used, via a simple matrix product to form an overall equivalent 
thermal transmission matrix for the entire fin. This thesis develops a thermal 
transmission matrix for the radiating rectangular, trapezoidal, and triangular fins both 
for the free space and non-free space environments. Test cases have been run and 
their solutions exactly match those contained in the literature. The thesis concludes 
with optimization studies for each profile considered where it is observed that simple 
algebraic equations can be employed to describe the optimum geometry. 
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I. INTRODUCTION 


Since the advent of the space age, extended surfaces have been extensively 
studied as a means for heat rejection from all sorts of vehicles. Because weight 
optimization is of the utmost importance in spacecraft design, it is desirable to have 
an accurate technique that will rapidly optimize the geometry of the extended surface 
used for rejecting the heat absorbed on or generated within the vehicle. To this point 
in time, the optimization problem involved the time consuming solutions of second 
order non-linear differential equations, by numerical methods or by cumbersome and 
possibly inaccurate graphical methods. A more advantageous method for determining 
the minimum mass, maximum radiant heat rejection extended surface geometry has 
been needed. 

The purpose of this study is to employ the new parameterization called the 
thermal transmission matrix to represent a single fin to evaluate the performance of 
several longitudinal fins. Each fin studied is subdivided into 100 subfins and the 
thermal transmission matrix is calculated for each subfin. Then, by a simple matrix 
multiplication of the individual subfin thermal transmission matrices, the overall 
transmission matrix is computed. Not only will this method produce the base and tip 
conditions of the overall fin entity, but the temperature and heat flow at any point 
between the fin base and fin tip can be determined. This iterative method considers 
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the spacecraft radiant heat rejection problem in both the free and non-free space 
environments. Optimization must be considered and simple algebraic equations that 
describe the fin geometry are formulated for the longitudinal fins of rectangular, 
trapezoidal, and triangular profiles. 

This new approach to fin synthesis and analysis permits an exact solution to any 
degree of accuracy that the user wishes to specify and produces the result in a 
computational efficient manner. 
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n. THEORETICAL PRESENTATION 


A. GENERALIZED DIFFERENTIAL EQUATION 

There are three distinct modes by which heat flows from one point to another. 

A rather simplistic definition of each can be stated. Conduction is the mode by which 
heat flows through a material by molecular motion. If the heat flows from a 
confining surface to a moving fluid, the mode is known as convection. If no 
intervening medium is present and the thermal energy is transported from one point to 
another in the form of electromagnetic waves, the mode is known as radiation or 
radiant heat transfer. 

The flow of heat by conduction is proportional to the temperature gradient and 
cross sectional area perpendicular to the direction of heat flow. 


^cd “ ■ 


dUx) 

dx 


( 2 . 1 ) 


Insertion of a proportionality constant yields 


A(x)^ ( 2 . 2 ) 

where the minus sign assures a positive heat flow in the presence of a temperature 
gradient which must be negative. This serves to define the thermal conductivity of 
the material written as k(T) because it is a material property almost always a function 
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k(T) = 


(2.3) 


^cd 


- A(x) 


dT{x) 

dx 


of temperature. Thus, the heat flow at x normal to a surface of area A(x) at a 
temperature T(x) is written as 


Qcd = ■ ^(^) 


dTlx) 

dx 


(2.4) 


In convection the heat flow from a confining surface at temperature T(x) to a 
fluid whose bulk temperature is will be proportional to the temperature difference 
and the amount of surface area. 


« S(x) [ Tlx) - T ] 


Insertion of a proportionality constant permits the representation 


= K ‘5(i^) t T(x) - ] 


and this serves to define the convective heat transfer coefficient h<,v as 


(2.5) 


( 2 . 6 ) 


^ _ _ ^_CV _ 

" ■ S(x) [ T(x) - ] 


(2.7) 
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and thus Newton’s Law of Cooling is written as 


9cv = Six) [ Tix) - 7 ; ] (2.8) 

where the convective heat transfer coefficient h,^ is written as a function of position 
and temperature. 

Radiation requires no intervening medium and, for a perfect emitter and 
absorber in full view of one another, the familiar Stefan-Boltzmann law shows that the 
heat flow is proportional to the surface area and the difference of the fourth powers of 
the absolute temperature. 

= a Six) [ - r/ ] (2.9) 

The proportionality constant used here is the Stefan-Boltzmann constant, 
a = 5.66961 x 10'* watt / m^ K'*. 

If the surfaces are not perfect emitters and not perfect absorbers and they are 
not in full view of one another, the departure from perfection is via the emissivity 
factor F, and the shape or arrangement factor F^. 

= o F^F^Six) [T/ - r/] (2.10) 

The problem of a finned surface radiating to outer space simplifies the foregoing 
equation (2.10). The temperature of the fm becomes T(x) and, because outer space is 
at about 2 K, T 2 can be assumed to be equal to 0. [Compare T(x)'’ = 200'* to T/ = 
2'*]. Moreover, F^ is equal to 1, the value of a small body in a large enclosure (small 
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Figure 2.1 Energy Balance 

fin in a large enclosure, outer space). Finally the emissivity factor for this case 
(small body in a large enclosure) is = €(x,T). The emissivity e is a property of the 
radiating fin surface and temperature. Thus equation (2.10) is rewritten as 

q^ = a e(x,T) S(x) T(xY (2-11) 

The difference between the rates of heat entering and leaving the element in 
Figure 2.1 by conduction (dq^j b - dq^j J must, in the steady state, be equal to the 
rates of heat dissipated from the face, L dx by convection (dq^vi and dq^v?) 
radiation (dq,] and dqr 2 ) plus any rates of heat that arrive from external sources (dq^ 
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and dqij). 

The presence of two terms for dq^^,, dq„ and dqj permit unequal values on the 
two faces. Thus with 


^cd.b 


dx 


k{x) b(x) L ^ 
dx 


dx 


( 2 . 12 ) 


the two convective dissipation terms are 


^cvi ^ ^cv 2 = Ldx[T 


+ h^,{x,T) Ldx[T-T,,] 


(2.13) 


and the two radiation dissipations terms with = 1 and F, = e(x,T) are 


^rl + = O ®(^»^ L dx [ ] 


+ o e(x,T) L dx 


4 >T ^ 


- 7 , 


52 


The energy balance with Tsi — Ts 2 — 0 is 


(2.14) 


_d 

dx 


Hx) 

dx 


^ = [ ^cv 2 (^’'^ ^cv 2 (^»^ ] L dx 7 


+ a e(x,T) L dx + q^j + q .2 (2-15) 

In space, if either face of the fin is subjected to an external heat input of the 


form 


q. = a L E dx 


(2.16) 
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where E consists of external heat inputs from solar, .Es, and/or terrestrial, Ep, sources 
then equation (2.15), the steady state thermal energy balance for a differential element 
dx, is rewritten as 


dx 


k(x) b{x) L 


dx 


dx = 




+ ^LT dx 


+ a e(x,T) LT^dx+laLEdx (217) 


B. LLMITING ASSUMPTIONS 

Further assumptions must be applied to the generalized differential equation 
(2.17) before a mathematical solution can be attempted. The assumptions employed 
by Murray' and Gardner^ are used for a starting point for further analysis. The 
Murray-Gardner Assumptions listed in Table 1 are essential in simplifying the 
formulation of the generalized differential equation (2.17) so that mathematical 
analysis is possible. Consider an arbitrary symmetrical longitudinal fm shown in 
Figure 2.2 of length (b-a), depth L, and width 6(x) dissipating radiant heat to a 
constant temperature T, surrounding environment. 

The profile area Ap is defined as 

Ap = J 6(x) dx (2.18) 

x=b-a 
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The initial conditions at the fin tip are 


T^ = nx=a) (2.19) 

and 

9a = 9( x=a ) (2.20) 

Applying the tenth Murray-Gardner assumption (Table 1), that no heat leaves from 
the tip of the fin, equation (2.20) leads to 

) . 0 ( 2 . 21 ) 

dx 

Two more boundary conditions are taken at the fin base 

nx=b) = T; (2.22) 

and 

9( x=b ) = (2.23) 

Applying the fifth and ninth Murray-Gardner assumptions (Table 1), allows one to 
ignore the radiant heat flow from the edges. The surface area S(x) perpendicular to 
the radiant heat flow is given by 

S(x) = 2 L Y (2-24) 

The face area S(x) through which the heat is dissipated depends on the arc length y, a 
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function of the fin length b. 


The cross sectional area A(x) perpendicular to the heat flow by conduction is 

Aix) = L b(x) (2.25) 



Because the second Murray-Gardner assumption (Table 1) states that the material is 
homogeneous, both the emissivity and thermal conductivity are assumed to be 
constant and uniform 
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TABLE 1 


THE MURRAY-GARDNER ASSUMPTIONS" 


1. The heat flow into the fin and the temperature at any point on the fin 
remains constant with time. 

2. The fin material is homogeneous, its thermal conductivity is the same in 
all directions and remains constant. 

3. The heat transfer coefficient between the fin and the surrounding medium 
is uniform and constant over the entire surface of the fin. 

4. The temperature of the medium surrounding the fin is uniform. 

5. The fin width is so small compared with its height that temperature 
gradients across the fin width may be neglected. 

6. The temperature at the base of the fin is uniform. 

7. There are no heat sources within the fin itself. 

8. Heat transfer to or from the fin is proportional to the temperature excess 
between the fin and the surrounding medium. 

9. There is no contact resistance between fins in the configuration or between 
the fin at the base of the configuration and the prime surface. 

10. The heat transferred through the outermost edge of the fin (the fin tip) 
is negligible compared to that through the lateral surfaces (faces) of the fin. 
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e{x,T) = e 


(2.26) 


k{T) = k (2-27) 

The contribution from convective dissipation from the faces of the fin is ignored 
because the arbitrary symmetrical longitudinal fin is operating in a space environment. 

C. GENERALIZED DIFFERENTIAL EQUATION FOR A FIN 

Combining the generalized differential equation (2.17) and equations (2.24) 
through (2.27), and assuming that there are no external heat inputs (qj, = q^ = 0), 
the steady state thermal energy balance for the differential element shown in Figure 
2.2 is 

—ik L - 2 L a e T{x)^ ^ = 0 (2.28) 

dx ) dx 

The temperature excess 0{x) is defined as the temperature difference between the 
temperature at point x and the constant temperature of the surrounding 
environment T,. 

e(x) = TXx) - (2.29) 

and from this, it is seen that 

dJXx) = dd{x) (2-30) 

Applying the temperature excess equations (2.29) and (2.30), the equation for the 
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steady state thermal energy balance for the differential element is rewritten as 



dx dx k dx 


(2.31) 


where h, is defined as the radiative heat transfer coefficient. Because = 0, 


(2.32) 


= a e 6(x)^ 


A fundamental concept of calculus is that the length of arc d 7 in a cartesian 
coordinate framework is 


(2.33) 


dy = [ dx^ + dy^] 


The derivative of the arc length is 



(2.34) 


The + or - sign is to be taken according as the arc length y increases or decreases as 
X increases. 

Combining the equation for the steady state thermal energy balance differential 
element (2.31) and the arc length equation (2.34), the generalized differential equation 
for a longitudinal fin of arbitrary symmetrical profile is 


^1 - lA e(:,) [l . '1^ = 0 (2-35) 

dx dx k 


Equation (2.35) is a second order, linear, homogeneous differential equation with the 
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boundary conditions 


and 


= 0 (2.36) 

dx 


0( x=b ) = 0^ (2.37) 


D. LONGITUDINAL FIN OF RECTANGULAR PROFILE 

The longitudinal fin of rectangular profile is illustrated in Figure 2.3. It should 
be noted that the origin of the x axis is located at the tip of the fin and not at the base 
of the fin. Because the width is constant over the entire profile, the width function 
6(x) is simply 

h(x) = 6, = (2.38) 

The profile area Ap is written as 

x=b 

= f 6(x) dx = 6i,b (2-39) 

1=0 


Inserting the width function (2.38) into the generalized differential equation for a fin 
equation (2.35) yields 


6 ,^ 
' dx^ 


k 


0(x) = 0 


(2.40) 
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or 


where 


- m2 0(x) = 0 (2.41) 

dx^ 


m 



(2.42) 


The radiative heat transfer coefficient h, is assumed to be a constant for a very small 
subfin where Tt » T^. 

The general solution of the generalized differential equation (2.41) is 


e(x) = Cj + C, (2.43) 

where C; and C 2 are arbitrary integration constants that are evaluated by applying the 
initial condition of equation (2.22). 


%ib) = 0^ = Cj e" ^ + C 2 e*'” * (2-44) 

The heat flow, qb, into the base of the fin is found by inserting the general 
solution of the generalized differential equation (2.43) and the cross sectional area 
perpendicular to the heat flow by conduction (2.25) into Fourier’s Law of 
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Figure 2.3 Longitudinal Fin of Rectangular Profile 


Conduction, equation (2.4). 

qix) = I it I ^[Cj e'"' + Cj (2-45) 

Simplifying and evaluating at x = b with the initial condition of equation (2.23), 
equation (2.45) becomes 

= L it /n ( C, e™ * - Cj e-'" * ) (2.46) 
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Solving for C, and C 2 yields 


Cl 


, ~m b 


0k + 


-m b 


2 K 




(2.47) 


and 


where 


b 


- 


0K - 


,m b 


2 y. 


<ib 


(2.48) 


= L [2 h^k b^) (2.49) 

is defined as the characteristic thermal admittance Yb of the fm. Likewise the 
characteristic thermal impedance, Zb, of the fm is the reciprocal of Yb 

^ ( 2 A, * (2-50) 

An examination of equations (2.43) through (2.48) reveals that only two of the 
four variables are independent and that the specification of any two variables allows 
for the determination of the remaining two. Combining equations (2.43) through 
(2.48) yields the two port inverse thermal transmission matrix for a longitudinal fm of 
rectangular profile 


■0. 


cosh mb 

-Z^ sinh mb 


% 



-Y^ sinh mb 

cosh mb 




(2.51) 
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where the inverse thermal transmission matrix is defined as 


T-* = 


cosh mb 
sinh mb 


-Zj sinh mb 
cosh mb 


(2.52) 


The thermal transmission matrix is defined as the inverse of the inverse thermal 


transmission matrix. 


T 


= = 


A 

C 


B 

D 


(2.53) 


where A, B, C, and D are the elements of the thermal transmission matrix. The 
thermal transmission matrix representation is 




cosh mb 

sinh mb 


X 



sinh mb 

cosh mb 




(2.54) 


E. LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE 

The longitudinal fin of trapezoidal profile is shown in Figure 2.4. Note that the 
origin of the axis is not at the tip of fin but located to the right of the tip 
where x = 0. 

The width function 5(x) for this profile is 

6(x) = 2 X tan ({) (2.55) 
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where </> is defined as 


4) = arctan 


6 . 


1 {b - a) 


(2.56) 


The profile area A is 


x-h 

Ap = f d(x) dx = ~ (6^ + 


(2.57) 


Inserting the width function (2.55) into the generalized differential equation for a fin 
[Equation (2.35)] yields 


or 


where 


2 X tan(<t)) 


d^Qjx) 

dx^ 


+ 2 tan((J)) 


dQ{x) 

dx 


0(4 = 0 (2.58) 

k cos(<})) 



dQjx) 

dx 


- m} 0(x) = 0 


(2.59) 


m 


1/2 


k sin(4)) 


(2.60) 


Equation (2.59) represents a second order, linear, homogenous differential equation. 
The general solution is expressed in terms of modified Bessel functions. The general 
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Figure 2.4 Longitudinal Fin of Trapezoidal Profile* 


solution to equation (2.59) is 

e(x) = Cj /j 2 w X ^ C^Kll m (2.61) 

The following transformation is made 

u = 2 m (2.62) 

such that at X = b 

u, = 2mb^f^ (2.63) 
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and at X = a 


= 2ma^^ (2.64) 

After applying the boundary conditions of equations (2.36) and (2.37) and 
making the transformations of equations (2.63) and (2.64), the general solutions 
become 

% = C, /„(u,) ♦ C, (2.65) 

and 

2 k b.L 

<it = -^- pi W - ^2 ^,(«^] (2.66) 


Simultaneous solution of the general solution equations (2.65) and (2.66) yield 


^1 “fe 


w % 


2 k by L 




(2.67) 


and 


^2 


w % - 




2 k bf^ L 




( 2 . 68 ) 


Substituting and rearranging the general solution equations (2.65) and (2.66) and the 
equations for C, and C 2 (2.67) and (2.68), the elements of the thermal transmission 
matrix for the longitudinal fin of trapezoidal profile are written as: 
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A =u^[ l,{uj KJiu^ ^ />^ K,(uJ] 


(2.69) 


B = 


2 k 


- 4 («.) W] ( 2 . 70 ) 


C = 2 k bL 

a 



[ /j(u^ K^iuJ - I^iUJ K^(u^] (2.71) 


D = u,l l,(uj K,iu^ + K,iuJ] 


(2.72) 


F. LONGITUDINAL FIN OF TRIANGULAR PROFILE 

The longitudinal fin of triangular profile shown in Figure 2.5 is a special case of 
the longitudinal fin of trapezoidal profile. There are many similarities between the 
longitudinal fin of triangular and trapezoidal profiles. The differential 
equation for the longitudinal fin of triangular profile is the same as the differential 
equation for the longitudinal fin of trapezoidal profile equation (2.59). 



dx^ dx 


(2.73) 
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Figure 2.5 Longitudinal Fin of Triangular Profile* 


where 


m 


k sm(4>) 


1/2 


(2.74) 


and 


4 ) 


arctan 


A 

2 b 


(2.75) 
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The profile area Ap is 


x=b 

= f 6(x) dx = ^ 6^ 


(2.76) 


However the evaluation of the arbitrary constants from the boundary conditions 
differ between the triangular and the trapezoidal profiles. As the value of x 
approaches zero the temperature excess at the tip 0^ must be equal to some finite 
value, therefore 

4( i/=0 ) = 1 (2.77) 

In order to physically obtain a finite temperature at x^ = u^ = 0, C 2 in equation 
(2.61) must be equal to zero because K2(0) is unbounded. The unbounded condition 
allows cancellation of the Ko(u) terms from the general solutions for temperature 
excess and heat flow. 

After cancellation and rearranging the general solutions for the longitudinal fin 
of trapezoidal profile equations (2.65) and (2.66), the general solution for the 
longitudinal fin of triangular profile becomes 


= 


2 k h^Lni^ 


W 


(2.78) 


As a result of the zero tip width, the edge parameters do not depend on the base 
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parameters. The thermal input admittance element for the longitudinal fin of 
triangular profile is 


Y = Si = 

' e. 


2 k /,(wp 


(2.79) 
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m. CASCADE ANALYSIS ALGORITHM 


A. GENERALIZED TWO PORT TRANSMISSION LINE MODEL 

The radiating fin is a distributed heat flow configuration which is analogous to 
the transmission line which is a distributed electrical network. The general two port 
transmission line model is shown in Figure 3.1. The model is described by two input 
terminals. At point "a" there is an applied voltage V, and a current I, flows into the 






Transmission 


'a V. 


Vb '‘» 


Line 






Figure 3.1 Two Port Transmission Line Model 
model. There is also a pair of terminals at point "b" where a voltage occurs and a 
current It flows out of the terminals. The two port model of the transmission line can 
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consist of any arrangement of passive elements within the two port but with no active 
current or voltage sources within the two port. 

It can be shown from network theory that the two port transmission line model 
in Figure 3.1 can be represented by the following equations: 

V^=AV^^BI^ (3.1) 


I^ = + D 

where the determinant of the ABCD matrix defined by 




A 

B 





C 

D_ 


It 


is equal to unity. 


A D - B C = 1 


(3.2) 


(3.3) 


(3.4) 


The values A, B, C, and D are constants and are defined as the elements of the 
transmission matrix. 

The electrical transmission line analogy can be directly applied to an equivalent 
thermal transmission line model where the heat flow is analogous to current and the 
temperature is analogous to voltage. Equations (3.1) and (3.2) are rewritten in the 
equivalent thermal transmission line model as: 

T^^AT^ + Bq^ (3.5) 
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= c + I? 


(3.6) 


B. CASCADE ALGORITHM 

The thermal transmission parameters are very useful in describing two port 
networks connected in a cascaded arrangement. It is an exercise in elementary two 
port theory. Examine the two equations for the two networks shown in Figure 3.2. 


and 


o 


A, 

Si' 




.^1 

1_ 





A, 

s; 




.^2 

1 _ 



(3.7) 


(3.8) 


Continuity for the configuration in Figure 3.2 dictates that 


and 


which can be represented as 


= 


( 3 . 10 ) 



(3.11) 


The overall thermal transmission parameter matrix for the two networks 
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cascaded in Figure 3.2 is by definition 


A 

B 




to 

B 2 

C 

D 


.^1 

^ 1 . 

.^2 

^ 2 . 



Figure 3.2 Cascaded Network 

and the thermal transmission matrix for a cascaded network of n networks is 

A B 
C D 

Equation (3.13) is an expression of what is referred to as the cascade algorithm. 

C. COMPUTER ANALYSIS 

The cascade algorithm expression (3.13) can be applied to a fin divided into n 
subfins each having the same length. Converting the temperature excess notation 



Ai B,' 


^2 ^2 


Bn 


q q 


^2 ^2, 
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back to the temperature notation (T^ « 0), the cascade algorithm expression (3.13) is 
written as 




A, B; 

■^2 ^2 







c, D, 


n 

9^. 


The base temperature Tb and the base heat flow are the conditions at the base and 
the tip temperature T, and tip heat flow q, are the conditions at the edge. 

A problem arises because the tip temperature T, and the base heat flow q^ are 
not known or not specified. The cascade algorithm expression (3.14) can not be 
solved directly because it is essentially a linear transformation that maps tip conditions 
to base conditions. However, a solution can be calculated indirectly. If one makes 
an initial estimation for the tip temperature T„ it can be used to calculate the radiant 
heat flow q, out of the tip of the fin. 

= o e 6^ L r/ (3.15) 

Thus T, and q, are the starting values for the computations using the cascade 
algorithm expression (3.14). At the conclusion of the process, a comparison 
involving all n subfins is then made between the known or specified base temperature 
Tb and the calculated base temperature. If the comparison indicates an error, the tip 
temperature T, is changed and the iterative process is repeated and continues until the 
error between the presumed known base temperature Tb and the calculated base 
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temperature are within an acceptable range. Figure 3.4 provides a complete flow 
chart of the cascade algorithm. 

Further amplification concerning the radiative heat transfer coefficient needs to 
be made. Recalling equation (2.32) where 

hj. = o e T(x)^ (3.16) 

it is not clear as to what temperature, T(x), should be applied to the computation of 
and hence to the cascade algorithm expression (3.14). The left and right faces of the 
differential element in Figure 3.3 have an average temperature. 

T + T. 

T = —2-* (3.17) 

■^av 2 



Figure 3.3 Differential Element 
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and as long as and Tt are not vastly different, equation (3.16) can be written using 


T 


h, = oeTj ( 3 . 18 ) 

The environmental parameter V is introduced for the case in which there are 
portions of the surrounding environment at some temperature other than T, » 0. An 
example of this is in the interception of thermal radiation from a nearby solar panel at 
some temperature not near absolute zero. The environmental parameter is defined as 


V = 


Ko 


2 o e 


( 3 . 19 ) 


where the parameter K 2 allows for the application of external heat inputs from solar, 
Es, or terrestrial, E^, sources or both. 


K^ = E a 


( 3 . 20 ) 


The environmental parameter is equal to zero for the free space case and to some 
number greater than zero but less than one for the non free space case. Substituting 
equation (3.19) into (3.18), the radiant heat transfer coefficient becomes 


hj. = a e 


^2 

2 r,. 


( 3 . 21 ) 


External heat inputs (qi and q^) from solar input Eg and terrestrial input Ej- are now 
applied to the cascade algorithm expression (3.14) through the application of the 
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environmental parameter V. 


D. PROGRAM DESCRIPTION 

The computer code consists of three analysis programs for each of the profiles 
considered, and a MAIN Program that is used for general housekeeping. The 
PERFORMANCE Analysis Program evaluates an entered set of parameters returning 
the characteristics (volume, profile area Ap, profile number PN, base heat flow q^, 
edge temperature T„ and the fin efficiency ri) of the fin along with the thermal 
transmission parameter matrix. The PROFILE AREA Analysis Program also 
evaluates an entered set of parameters but it returns a data set of 100 fins having a 
constant profile area Ap with the base fin width d and length b which are used as the 
input variables. The OPTIMIZE Analysis Program is similar to the 
PERFORMANCE Analysis Program with the exception that it returns the optimum 
fin geometry which results in the minimum mass and maximum base heat flow fin. 

The MAIN Program provides menu driven access to the various analysis 
programs. The user has the ability to set the computational main program error 
coefficient. The user also has the ability to set the save data switch which controls 
data flow into an ASCII formatted disk file. Commas are used for the data delimiter 
which allows file importation into MATLAB, Word Perfect, Draw Perfect, and Math 
CAD for additional analysis or for plotting. An additional feature provided is the unit 
conversion module. Table 2 list the available units that are provided to the user by 
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the MAIN Program. Figure 3.4 provides a flow chart of the cascade algorithm 
expression (3.14) which is common to all of the analysis programs. 

The PERFORMANCE Analysis Program requires the following input conditions 
to be entered: 

• base temperature 

• density 

• emissivity 

• environmental parameter 

• fin base width 

• fin edge width* 

• fin depth 

• fin length 

• thermal conductivity 

The program calculates and returns the following values: 

• base heat flow 

• edge temperature 

• fin efficiency 

• profile area 

• profile number 

• thermal transmission parameter matrix 

• total mass 

The fin efficiency tj and the profile number PN are extensively referred to in 
literature and are provided for use in making comparisons. The fin efficiency and the 
profile number are defined as 


2 a e jb L 


(3.22) 


'Applies for only the longitudinal fin trapezoidal profile. 
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Figure 3.4 Cascade Algorithm Flow Chart 
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TABLE 2 


UNITS AVAILABLE IN THE MAIN PROGRAM 


Admittance 

• watts / °C, BTU / hr °F 

Area 

• m^, cm^, fF, in^, mm^ 

Heat 

• watts, kw, BTU / sec, BTU / hr 

Impedance 

• °C / watts, °F hr / BTU 

Length and Depth 

• m, cm, ft, in, mm 

Mass 

• kg, grams, Ibm, oz 

Temperature 

• K, °C, °F, °R 

Thermal Conductivity 

• watts / m K, BTU / ft hr °F 

Width 

• m, cm, ft, in, mm 


PN 


2 a e 

F6^ 


(trapezoidal) 


(3.23) 


PN = 




o e b^ (rectangular) (3.24) 
k 6 ^ 


If the save switch is switched on, the following differential element conditions are 


48 







written to a disk file: 


• base heat flow 

• base temperature 

• distance from the tip of the fin 

The OPTIMIZE Analysis Program requires the following input conditions to be 
entered: 


• base temperature 

• density 

• emissivity 

• environmental parameter 

• fin edge width^ 

• fin width 

• thermal conductivity 

The program calculates and returns the following optimum values: 

• base heat flow 

• edge temperature 

• fin base width 

• fin efficiency 

• fin length 

• profile area 

• profile number 

• thermal transmission parameter matrix 

• total mass 

The program saves the same quantities that the PERPORMANCE Analysis Program 
saves depending on the condition of the save switch. 

The PROFILE AREA Analysis Program requires the following input conditions 


^Applies for only the longitudinal fin trapezoidal profile. 
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to be entered: 


• base temperature 

• density 

• emissivity 

• environmental parameter 

• profile area 

• starting fin base width value 

• ending fin base width value 

• taper ratio^ 

• fin depth 

• thermal conductivity 

The program calculates and writes to a disk file the following parameters for 100 fins 
having the same profile area Ap within the specified width range: 

• base heat flow 

• base width 

• base temperature 

• edge temperature 

• fin efficiency 

• thermal transmission parameter matrix 

• total mass 

The analysis programs all use the thermal transmission matrices computed in 
Chapter 2 in the cascade algorithm with the exception of the longitudinal fin of 
triangular profile. The triangular and trapezoidzil profile analysis are very similar. 

The analysis program uses the thermal admittance equation (2.79) to compute the base 
conditions for the first subfin having a thickness of 5.. The sub fin is located at the tip 
of the fm. The base conditions are then applied and used as the tip conditions in the 
thermal transmission matrix for the trapezoidal profile (2.69) through (2.72). The 


^Applies for only longitudinal fin of trapezoidal profile. 
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cascade algorithm is then used to calculate the thermal transmission matrices for the 
remaining 99 subfins in order to finally determine the overall thermal transmission 
matrix and the base conditions. The trapezoid that is being used in the calculation has 
an almost zero tip width compared to the base width such that 

Sfc = 100 6^ (3.25) 

Borland Turbo C4-+ was used to write the computer code. The computer code 
listing appears in Appendix A. Full compiler optimization was selected in order to 
decrease the run times. Average run times on an IBM 486/487 machine appear in 
Table 3. 


TABLE 3 

AVERAGE PROGRAM RUN TIMES 



PERFORMANCE 

(sec) 

PROFILE AREA 
(sec) 

OPTIMIZE 

(sec) 

Rectangular 

4 

97 

4 

Trapezoidal 

4 

600 

4 

Triangular 

5 

625 

5 
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IV. PROGRAM ANALYSIS RESULTS AND DISCUSSION 


A. PROGRAM ERROR TOLERANCE 

In order to allow for meaningful interpretation, an understanding of the accuracy 
of the results produced from the cascade algorithm is vital. A test to show 
convergence is used to illustrate the approximate relative error. 

Table 4 provides an illustration of how the error coefficient in the MAIN 
Program affects the accuracy of the program results. Table 4 is produced by 

TABLE 4 

PROGRAM ERROR TOLERANCE 


Program 


Error 

Ta 

« % 

qb 

« % 

Coeff % 

(K) 

error 

(watts) 

error 

1.000 

282.18 

0.0106 

926.28 

2.0918 

0.500 

282.24 

0.0106 

931.83 

1.5052 

0.100 

282.19 

0.0071 

946.36 

0.0307 

0.050 

282.22 

0.0035 

946.94 

0.0920 

0.010 

282.22 

0.0035 

946.10 

0.0032 

0.005 

282.21 

0.0000 

946.05 

0.0021 

0.001 

282.21 

0.0000 

946.07 

0.0000 


running the same input parameters many times while decreasing the error coefficient 
in the MAIN Program. The PERFORMANCE Analysis Program for the longitudinal 
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fin of rectangular profile was used. The Table 4 results are applicable for any of the 
analysis programs because of the commonality of the cascade algorithm. As the error 
coefficient in the MAIN Program decreases, the output value converges. The 
approximate relative error e, is defined as 




True Value - Approximate Value 
True Value 


100 


(4.1) 


The exact true value is unknown, however the value produced when the error 
coefficient in the MAIN Program is equal to 0.001 can be assumed to be 
approximately equal to the exact true value. 

Examining the results illustrated in Table 4, one can conclude that when the 
error coefficient in the MAIN Program is less than or equal to 0.1 the approximate 
relative error observed in the output is generally much less than the error coefficient 
in the main program. 


B. GRAPHICAL EXAMPLES 

The purpose of this section is to describe and provide examples on how each of 
the three analysis programs can be used. The computer code listings are provided in 
Appendix A. 

1. Example One: Longitudinal Fin of Rectangular Profile Analysis 

A longitudinal fin of rectangular profile is 0.5 meter long, 1.0 centimeter 
thick, and 1.0 meter deep. The fin is made of aluminum (k = 209.4 watt / m K, 
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p = 2.6 kg / dm\ and c = 0.85). If the base temperature is 400 K, what are the tip, 
midpoint, and base conditions (error < 0.1%) both in the free space and non free 
space environments? 

Because the geometry of the fin is already known, the PERFORMANCE 
Analysis Program is used for the evaluation. The save switch in the MAIN Program 
must be placed in the save position in order to later retrieve the midpoint conditions 
from the disk file. The computer output display for the free space environment 
(V=0) appears in Table 5. 

Figures 4.1 and 4.2 provide a detailed graphical representation of the 
required conditions both in the free space and non free space environments. The base 
conditions are read from where the length is equal to 50 centimeters and the tip 
conditions are read from where the length is equal to 0 centimeter. In the free space 
environment the tip, midpoint, and base conditions for this particular fin geometry are 
listed in Table 6. 
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TABLE 5 


PERFORMANCE ANALYSIS OUTPUT DISPLAY 


THE LONGITUDINAL FIN 

OF 

RECTANGULAR PROFILE 

Density 

= 

2.60 kg / dm^3 

Fin Depth 

= 

1.00 m 

Total Mass 

= 

13.00 kg 

Fin Length 

= 

0.50 m 

Fin Width 

= 

1.00 cm 

Emissivity 

= 

0.85 

Profile Area 

= 

50.00 cm^2 

Profile Number 

= 

0.607 

Base Heat Flow 

= 

732.85 watts 

Fin Efficiency 

= 

0.594 

Base Temperature 

= 

400.07 K 

Edge Temperature 

= 

324.45 K 

Thermal Conductivity 

= 

209.40 watts / m K 

Environmental Parameter 

= 

0.00 

Transmission Parameter Matrix 

A = 1,2287 

B 

1 = 0.2583 deg C / watt 

B = 2.2374 watt / deg C 

D 

1 = 1.2842 


TABLE 6 

FIN CONDITIONS IN FREE SPACE 


Location 

Temp 

q 


(K) 

(watts) 

Base 

400.07 

732.85 

Midpoint 

341.66 

292.33 

Tip 

324.45 

10.68 
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Figure 4.1 Heat Dissipation From A Rectangular Profile Fin 
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Figure 4.2 Temperature Distribution For A Rectangular Profile Fin 
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Equation (3.4) is used below to verify that the determinant of thermal 


transmission matrix listed in Table 5 is indeed equal to unity. 


det 


1.2287 0.2583 

2.2374 1.2842 


1.0000 


(4.2) 


2. Example Two: Longitudinal Fin of Rectangular Profile Area Analysis 
A longitudinal fin of rectangular profile 1.0 meter deep is made of 
aluminum (k = 209.4 watt / m K, p = 2.6 kg / dm^, and e = 0.85). It is required 
to remove a minimum of 500 watts of heat from the base of a fin at a temperature of 
400 K. Determine the geometry of the fin using several values for the profile area 
Ap- 


The PROFILE AREA Analysis Program is used because the geometry of 
the fin is unknown. Figure 4.3 provides a graphical representation which shows the 
constant profile area Ap relationship as a function of width versus heat flow. Any of 
the dimensions listed in Table 7 would meet the minimum requirement to remove the 
500 watts of heat from the base of the fin. 

Examining Figure 4.3, one can see that the maximum heat flow occurs at 
a particular width value. More heat is dissipated as the profile area Ap is increased. 
The value of the width at the point of maximum heat flow also increases as 
the profile area Ap increases. The point of maximum heat flow is considered as the 
optimum geometry (maximum heat flow, minimum mass) for the particular profile 
area Ap. A summary of the optimum geometry points extrapolated from the data files 
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FIGURE 4.3 Profile Area Analysis 
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TABLE 7 


PROHLE AREA ANALYSIS POINTS 


Ap 

5b 

b 

Qb 

(cm^) 

(cm) 

(m) 

(watts) 

25 

0.334 

0.748 

507.4 

100 

0.490 

2.041 

633.7 

200 

0.724 

2.762 

771.2 

300 

0.334 

8.982 

523.1 


produced from the PROFILE AREA Analysis Program appear in Table 8. 

TABLE 8 

OPTIMUM GEOMETRY POINTS FOR THE RECTANGULAR PROFILE 


Ap 

5bo 

bo 

MasSo 

Qbo 

qbo/m 

(cm^) 

(cm) 

(m) 

(kg) 

(watts) 

(w/kg) 

25 

0.646 

0.387 

6.50 

579.76 

89.19 

100 

1.543 

0.648 

26.00 

924.53 

35.56 

200 

2.440 

0.819 

52.00 

1165.83 

22.42 

300 

3.103 

0.967 

80.00 

1334.40 

16.68 


3. Example Three: Comparison of Various Profiles 

A longitudinal fin 1.0 meter deep having a profile area Ap equal to 50 cm^ 
is made of aluminum (k = 209.4 watt / m K, p = 2.6 kg / dm^ and 
e = 0.85). The base temperature is held at 400 K. Compare the relative 
performance between the longitudinal fins of rectangular, trapezoidal, and triangular profiles. 
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The same approach is used as in the previous example. A longitudinal fin 
of trapezoidal profile having a taper ratio TR = 0.5 is used. The taper ratio is 
defined as 

TR = — (4.3) 

Figure 4.4 provides a graphical representation which shows the relationship between 
the three profiles being considered. 

The curve in Figure 4.4 produced by the longitudinal fin of trapezoidal 
profile is midway between the curves produced by the longitudinal fins of rectangular 
and triangular profiles. Of more importance is the fact that at the location of the 
optimum geometry points, the longitudinal fin of triangular profile produces the 
largest heat flow while the longitudinal fin of rectangular profile produces the smallest 
heat flow. The optimum geometry points extrapolated from the disk file generated 
from the PROFILE AREA Analysis Program appear in Table 9. Not only is there an 
optimum fin geometry for a particular profile area Ap but there is also an optimum 
profile shape which produces maximum heat flow. 

4. Example Four; Optimum Geometry Comparison 

A longitudinal fin 2.0 meters deep (e = 0.95) is required to remove at 
least 1000 watts from a spacecraft in a free space environment (V=0). The base 
temperature is held at 400 K. Using the list of materials appearing in Table 10 
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Figure 4.4 Comparisons of Various Profiles 
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TABLE 9 


OPTIMUM GEOMETRY POINTS FOR THE VARIOUS PROFILES 


Profile 

RECT 

TRAP 

TRI 


Taper 

Ratio 

1.0 

0.5 

0.0 


5 . 

(cm) 

0.958 

1.270 

1.309 


b 

(m) 

0.522 

0.524 

0.763 


732.47 

772.93 

826.78 


Qb/m 

(W/kg) 

56.343 

59.456 

63.598 


Mass Qb 
(kg) (watts) 

13.00 
13.00 
13.00 


conduct a heat flow per unit mass comparison using the longitudinal fins of 
rectangular and triangular profiles. 


TABLE 10 

MATERIALS USED IN THE OPTIMUM COMPARISON* 



Chemical 

P 

k 

Material 

Symbol 

kg/dm^ 

w / m K 

Aluminum 

A1 

2.6 

209.4 

Copper 

Cu 

8.8 

384.0 

Iron 

Fe 

7.2 

58.0 

Magnesium 

Mg 

1.7 

157.0 

Molybdenum 

Mo 

10.2 

145.0 

Silicon 

Si 

2.3 

83.0 

Silver 

Ag 

10.5 

407.0 

Zinc 

Zn 

6.8 

140.0 
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The OPTIMIZE Analysis Program is used for the analysis. The 
parameters were entered individually into the OPTIMIZE Analysis Program and the 
results are listed in Table 11. 


TABLE 11 

OPTIMUM COMPARISON RESULTS 




«bo 


MasSg 

qbo 

qbo/ni 

Material 

cm 

cm^ 

kg 

watts 

w / kg 

Ag 

RECT 

0.40 

21.73 

45.63 

997.86 

21.86 

Ag 

TRI 

0.54 

15.09 

31.69 

995.45 

31.41 

A1 

RECT 

0.79 

42.24 

21.96 

998.96 

45.49 

A1 

TRI 

1.05 

29.33 

15.25 

995.47 

65.27 

CU 

RECT 

0.43 

23.03 

40.54 

997.94 

24.61 

Cu 

TRI 

0.57 

15.99 

28.15 

995.45 

35.36 

Fe 

RECT 

2.84 

152.49 

221.11 

1003.99 

4.54 

Fe 

TRI 

3.79 

105.89 

153.54 

995.75 

6.48 

Mg 

RECT 

1.05 

56.33 

19.60 

999.74 

51.01 

Mg 

TRI 

1.40 

39.12 

13.61 

995.49 

73.14 

Mo 

RECT 

1.14 

61.00 

124.43 

1000.00 

8.04 

Mo 

TRI 

1.51 

42.36 

86.41 

995.54 

11.52 

Si 

RECT 

1.99 

106.56 

49.66 

998.41 

20.11 

Si 

TRI 

2.65 

73.99 

34.48 

995.59 

28.87 

Zn 

RECT 

1.18 

63.17 

85.92 

1000.12 

11.64 

Zn 

TRI 

1.57 

43.87 

59.66 

995.50 

16.68 
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Figure 4.5 shows a bar graph representation of the heat flow per unit mass 
relationship using the various materials listed in Table 10. It is evident from Figure 
4.5 that the optimized longitudinal fm of triangular profile made of magnesium is the 
most efficient material considered of those listed in Table 10, with aluminum being a 
close second. The melting temperature of both materials (= 900 K) is far above the 
base temperature (400 K) stated in the problem statement. Comparing Tables 10 and 
11 materials having a small value for the density (p < 3.0 kg / dm^) and a large 
value for the thermal conductivity ( k > 150 watt / m K) produce large values for 
heat flow per unit mass. 

C. DERIVATION OF WORKING OPTIMUM GEOMETRY EQUATIONS 

As shown in Tables 8 and 9 optimum geometry points can be determined using 
the PROFILE AREA Analysis Program. The overall success of finding the 
dimensions for the optimum geometry using the PROFILE AREA Analysis Program 
has several limitations. The first limitation is that the starting and ending fin base 
width values must be known. And it can be argued that a second limitation is that 
excessive analysis program run times which might exceed ten minutes. 

Relationships for the optimum geometry for the longitudinal fins of rectangular, 
trapezoidal, and triangular profiles would remove both limitations. Kern and Kraus’ 
and Chung and Nguyen'® stated that the dimensions for optimum fin geometry can 
be expressed by 
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Figure 4.5 Optimum Profile Comparison Using Various Materials 
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(4.6) 


where Pi(TR), P 2 (TR), and P 3 (TR) are defined as polynomials in which the taper ratio 
TR is the independent variable. 

In general the polynomial P is expressed in the form of 


PiTR) = Cj TR TR + ... + C„ TR ° (4.7) 

where n is the degree of the polynomial. 

The polynomials Pj, P 2 , and P 3 are determined from running the PROFILE 
AREA Analysis Program using a slightly modified code. The computer code listing 
appears in Appendix B. A random number generator is used to generate the input 
parameters. The random number generator produces input parameters in the range 
listed in Table 12. Murray-Gardner (Table 1) assumption number five (the fin width 
is small compared with its depth and length) is checked in order to insure that the 
original assumptions are kept valid. 
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TABLE 12 


RANGE OF THE RANDOM INPUT PARAMETERS 


0.01 < e < 0.99 
20 < L ( cm ) < 220 
10 < k ( watt / m K) < 210 
100 < Tb ( K ) < 600 
25 ^ Qb ( watts ) < 1025 
0.01 < Taper Ratio TR < 0.99 


The program enters into a loop calculating the heat flow qb using the random 
input parameters. The program remains in this loop until the maximum heat flow Qb 
is found. Once found the program writes the random input parameters and the value 
determined to be the maximum heat flow qb to a disk file. Another set of random 
input parameters is generated and the maximum heat flow qb once again is calculated. 
This sequence continues until 100 optimum data sets are written to a disk file. 

Three data sets using the random input parameters are produced, one each for the 
longitudinal fins of rectangular, trapezoidal, and triangular profiles. 

The last step in determining the polynomials Pi, Pj, and P 3 is to fit the collected 
data stored in the disk file to a "best" polynomial fit having the smallest degree d 
using the least square principle. Equations (4.4) through (4.6) are rewritten as 
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(4.8) 


P(TR) 


Lb^oe T; 

<lb 



(4.9) 


PiTR)^ = APo k t; P 


(4.10) 


Equations (4.8) through (4.10) are placed into a MATLAB "m" file. Each of the 
collected data files are read into MATLAB and were evaluated using the MATLAB 
polynomial fitting function polyfit with the taper ratio TR as the independent variable 
using varying degrees d in order to find the polynomials P,, P 2 , and P 3 which "best" 
approximate the previously collected data from the disk file. The results appear in 
Table 13. The listing of the MATLAB evaluation m files appear in Appendix C. 


The standard error e, is defined as 


1/2 


E (>’. - y) 


(4.11) 


N - 1 


where y; is the computed value of y; and there are N values. 
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TABLE 13 


COEFFICIENTS OF THE OPTIMUM POLYNOMIALS 



Degree 




Standard 

Profile 

d 

c, 

C: 

C3 

Error 

PI RECT 

2 

0.8675 

0.0000 

0.0000 

3.54e-3 

PI TRAP 

2 

-0.0290 

0.8138 

0.0000 

2.12e-3 

PI TRI 

2 

0.0000 

0.9040 

0.0000 

2.55e-3 

P2 RECT 

2 

1.8648 

0.0000 

0.0000 

1.63e-2 

P2 TRAP 

3 

0.3881 

-1.0324 

2.7522 

1 .22e-2 

P2 TRI 

2 

0.0000 

2.4852 

0.0000 

1.40e-2 

P3 RECT 

2 

1.6178 

0.0000 

0.0000 

2.06e-2 

P3 TRAP 

2 

0.5147 

1.1382 

0.0000 

1.24e-2 

P3 TRI 

2 

0.0000 

1.1234 

0.0000 

9.46e-3 


The equation for the optimum profile area Ap,, for either the rectangular, 
triangular, or trapezoidal profiles is written as 


= ( 1.1382 + 0.5147 TR ) 


/ \3 

1 




L T\ 



(4.12) 
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V. CONCLUSIONS 


This thesis has described a new technique that uses a parameterization called the 
thermal transmission matrix to represent a single fin to evaluate the performance of 
several fins. The algorithm is based on the simple ideal of subdividing the fin into 
100 subfins and then the thermal transmission matrix is calculated for each subfin. 
Then, by a simple matrix multiplication of the individual subfin thermal transmission 
matrices, the overall thermal transmission matrix is computed. Not only will this 
produce the base and tip conditions of the overall fin entity, but the temperature 
distribution and heat flow at any point between the fin base and the fin tip can be 
easily determined. 

Several conclusions can be made from the Chapter 4 Program Analysis Results 
cuid Discussion. It was observed that the profile area Ap, increases as the cube of the 
heat flow q,,. In order to double the heat flow q,,, the profile area Ap must be 
increased eight times. 

Material selection is very important. The weight of the fin is proportional to 
the density and the thermal conductivity is inversely proportional to the profile area. 
Materials having a thermal conductivity greater than 150 watt / m K and a density 
less than 3 kg / dm^ are generally the best suited materials. For the same material, 
environmental conditions, heat flow, and base temperature, the optimum longitudinal 


71 


fin of triangular profile requires only 69 percent as much material as the optimum 
longitudinal fin of rectangular profile, with the optimum longitudinal fin of 
trapezoidal profile falling in between. 

Simple working relationships to compute the optimum fin dimensions are found. 
Table 16 compares the findings from this thesis with the results from previous work 
by others. The differences between the optimum dimension polynomial coefficients 
listed in Table 16 are very small, which verifies the accuracy of the thermal 
transmission matrix parameterization technique. 

TABLE 16 

OPTIMUM DIMENSION POLYNOMIAL COEFFICIENTS 


RECTANGULAR 


TRIANGULAR 


PI P2 P3 


.8675 

1.8648 

1.6178 

This Thesis 

.8846 

1.8473 

1.6341 

Chung and Nguyen’ 

.8842 

1.8579 

1.6056 

Liu’^ 

.8846 

1.8473 

1.6341 

Bartas and Sellers‘S 

.9040 

2.4852 

1.1234 

This Thesis 

.9485 

2.3104 

1.0957 

Chung and Nguyen 

.9546 

2.3027 

1.0991 

Nilson and Curry^^ 

.9091 

2.4255 

1.1025 

Kern and Kraus’^ 


This new approach to fin synthesis and analysis permits an exact solution to any 
degree of accuracy that the user wishes to specify and produces the result in a 
computational efficient manner. 
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AI’PENDIX A COMPUTER CODE LISTING 




Main Program 

A*************************************************** ******** *★***«***★★ -k j 


#include <process.h> 
#include <conio.h> 
^include <stdlib.h> 
^include <stdio.h> 

^define Yes 1 
^define No 0 


char 

*S1 C] 

= 

L 

M 

watt / deg C", " BTU / hr deg F" }; 


char 

*S2C] 

= 

L 

M 

m " 2" 

, " cm " 2", " ft " 2", " in " 

2", '* m\ " 2" } 

char 

*S3[] 

= 

L 

II 

watts" 

, " kw", " BTU / sec", " BTU / 

hr" >; 

char 

*S4[] 

= 

L 

II 

deg C 

/ watt", " deg F hr / BTU" ); 


char 

*S5[] 

= 

c 

II 

m", " 

cm", " ft", " in", " mm" }; 


char 

*S6[] 

= 

c 

II 

kg", " 

grams", " Ibm", " oz" }; 


char 

*S7n 

= 

i 

II 

deg K" 

, " deg C", " deg F", " deg R" 

}; 

char 

*S8[] 

= 

i 

II 

watt / 

m deg K", " BTU / ft hr deg F" 

' ); 

char 

*s9n 

= 

i 

II 

No", " 

Yes"}; 


char 

*0[] 

= 

i 


■‘<1> 

watt / deg C <2> BTU / hr deg 



"<1> m " 2 <2> cm " 2 <3> ft ^ 2 <4> in ^ 2 <5> mm ^ 2'*, 

'*<1> watts <2> kw <3> BTU / sec <4> BTU / hr", 

'*<1> deg C / watt <2> deg F hr / BTU", 

"<1> m <2> cm <3> ft <4> in <5> mm", 

"<1> kg <2> grams <3> Ibm <4> oz", 

"<1> deg K <2> deg C <3> deg F <4> deg R", 

"<1> watt / m deg K <2> BTU / ft hr deg F", 

"<1> No <2> Yes" }; 

void print_main_menu(), print_menu_one(), print_menu_two(); 

void setupO, overlaysO, update_units(); 

float error=0.1; 

int T=0, L=0, M=0, S=0, H=0, d=4, K=0, A=1, 2=0, Y=0; 

char main_key=0; 

int main()< 

FILE *data_in; 

if ( (data_in = f open("uni ts.ext", "r")) != NULDL 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&A,&2,&Y,&error); 
fclose(data_in); } 

print_main_menu(); 

doL 

V if (main_key == '1*) setupO; 

if (main_key >= '2' 8,8, main^key <= '4') overlaysO; } 

while (( main_key = getchO ) != '5'); 

windowC1,1,80,25); 

clrscr(); 

_setcursortype(_NORMALCURSOR); 
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return 0;> 


void over I ays(void){ 
char key=0; 

int Quit = No, result; 
prin t_menu_one(); 
do 

switch (key){ 


case 'A': clrscrO; 
if (main_key 
if (main_key 
if (main_key 
break; 


Quit = Yes; 
s= 12*) result 
== '3*) result 
== '4‘) result 


spawnl(P_WAIT, "RECT.EXE", NULL); 
spawnl(P_WAIT, "TRAP.EXE", NULL); 
spawnl(P_WAIT, "TRI.EXE", NULL); 


case 'B': clrscrO; Quit = Yes; 


if 

(main_key == ‘2') 

result 

if 

(main_key == *3') 

result 

if 

(main key == '4') 

result 

break; 



spawnl(P_WAIT, "ARECT.EXE", NULL); 
spawnl(P_WAIT, "ATRAP.EXE", NULL); 
spawnl(P_WAIT, ••ATRI.EXE", NULL); 


case 'C: clrscrO; 
if (main_key 
if (main_key 
if (main_key 


Quit = Yes; 

== '2*) result 
r= '3') result 
== '4') result 


spawnl(P_WAIT, "ORECT.EXE", NULL); 
spawnl(P_WAIT, "OTRAP.EXE", NULL); 
spawnl(P_WAIT, "OTRKEXE", NULL);} 


while (( key = toupper(getch()) ) != 'D' && Quit == No); 


print_main_menu() ; 
return;} 

void setup(void){ 
char key=0; 


print_menu_two(); 
do 

switch (key){ 


case 'A': 


case 'B': 


case 'C: 


case 'D': 


case 'E': 


case 'F‘: 


case 'G': 


gotoxy(2,17); clreolO; cprintf ("%s*',0[0] ); 
while (key != '1' && key != '2*) key = getchO; 

Y = key - '1';print_menu_two(); break; 

gotoxy( 1,17); clreolO; cprintf("%s",0[1] ); 

while (key != O' && key != *2'&& key != '3*&& key != '4'&& key != '5') 
key = getchO; 

A = key - '1 *;print_menu_two(); break; 
gotoxyd , 17); clreolO; cprintf ("%s",0[2]); 

while (key != '1' && key != '2'&& key != '3'&& key != *4') key = getchO; 
H = key - '1';print_menu_two(); break; 

gotoxy(2,17); clreolO; cprintf("%s",0[3] ); 
while (key != *1' && key != '2') key = getchO; 

Z = key - M';print_menu_two(); break; 

gotoxy(2,17); clreolO; cprintf("%s",0[4]); 

while (key != '1' && key != '2'&& key != •3'&& key != •4'&& key != *5') 
key = getchO; 

L = key - '1 *;print_menu_two(); break; 
gotoxy(3,17); clreolO; cprintf("7oS'\0[5) ); 

while (key 1= '1* && key != ■2'&& key != '3'&& key 1= '4') key = getchO; 
M = key - '1';print_menu_two(); break; 

gotoxyd , 17); clreolO; cprintf (“%s",0[6] ); 

while (key != '1' && key != '2'&& key != '3'&& key != '4') key = getchO; 


case 'H': 


case 'r: 


case 'J': 


case 'K*: 


T = key - '1 *;print_menu_two(); break; 

gotoxyC 1,17); clreolO; cprinrf(*'%s",0[7] ); 
while (key != '1‘ && key != '2') key = getchO; 

K = key - '1';print_menu_two(); break; 

gotoxy(2,17); clreoU); cprintf(‘'%s",0[4] ); 

while (key != 'I* && key != 'Z'&S key != key != •A'&S. key != 'S’) 

key = getchO; 

d = key - '1 *;print_menu_two(); break; 

gotoxy( 12,17); clreolO; cprintf("%s",0[8] ); 
while (key != '1' && key != ’2') key = getchO; 

S = key - '1 *;print_menu_two(); break; 

gotoxy(22,12); clreolO; scanf("%f",&error); 
print_menu_two();> 


while (( key = toupper(getch()) ) 1= 'L'); 


update^units(); 
pr1nt_main_menu(); 
return;} 


void update_units(void){ 

FILE *unit; 

unit = fopenC'units.ext","w"); 

fprintf(unit,"%d %d 7od %d %d %d %d %d %d %d %f",T,L,M,S,H,d,K,A,Z,Y,error); 

fclose(unit); 

return;} 

void print_main_menu(void){ 

_setcursortype(_NOCURSOR); 
windowd, 1,80,25); 
clrscr(); 

gotoxy(27,3); cprintf("EXTENDED SURFACE ANALYSIS"); 
window(20,6,80,25); 


cprintf("<1> 
cprintf (‘'<2> 
cprintf("<3> 
cprintf("<4> 
cprintf("<5> 
return;} 


Setup / Initialization\r\n"); 

Fin of Rectangular Profile\r\n"); 
Fin of Trapezoidal Profile\r\n"); 
Fin of Triangular Profile\r\n"); 
Return to DOS"); 


void print_menu_one(void){ 

clrscr(); 
gotoxy(1,1); 


cprintf (".\r\n"); 

cprintfC <A> Predict the Performance\r\n"); 

cprintf(" <B> Analyze the Profile Area\r\n"); 

cprintf(" <C> Optimize the Dimensions\r\n"); 

cprintf(" <D> Return to Main Menu\r\n"); 

cprintfC."); 

return;} 


void print_menu_two(void)< 


cl rscrO; 
gotoxy(1,1); 
cprintfC — 
cprintf("<A> 
cprintf C<B> 
cprintf("<C> 
cprintf("<D> 
cprintf C<E> 


Admittance 

Area 

Heat 

Impedance 

Depth 


: %s\r\n", S1[Y]); 
: %s\r\n", S2[A]); 
: %s\r\n", S3[H]); 
: %s\r\n", S4[Z]); 
: %s\r\n", S5[L]); 


\r\n"); 
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cprintf("<F> 
cprintf("<G> 
cpn'ntf ("<H> 
cprintf("<I> 
cprintf("<J> 
cprintf("<K> 
cprintf("<L> 
cprintf("— 
return;} 


Mass : %s\r\n'', S6[M]); 

Temperature : %s\r\n", S7[T]); 

Thermal K ; %s\r\n", SSTK]); 

Width : %s\r\n», S5[d]); 

Save Data : %s\r\n", S9[S]); 

Percent Error: %f\r\n", error); 
Return to Main Menu\r\n"); 
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j-k ★***★*★*★★★**★**•***★*★*********★★*****★*★**★*****★*★**★★*****★★★**********★ 


PERFORMANCE Analysis Program 
Longitudinal Fin of Rectangular Profile 

■k-kk-k-k-kkk-k'kkkkkk'k-kk-k-k'k-k-kk-kk-kkkk'kkkkkkk-k-k-k'k-kkkk'kk-k-k-k-k-k-k-k-k-k-k-k-kk-k-kit-k-kkk-k-k-kkiftck-k-k k j 


#incLude <stdio.h> 
^include <stdlib.h> 
#include <math.h> 

#include <conio,h> 
#include "units.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 100 

sigma 5.66961e-8 


char f i le_name[20]; 

double A_(), B_(), C_(), D_(), m(), hr(), Yb(), A, B, C, D; 

double length, width, depth; 

double inc, overall_diff = 100, Tmax = 0, Tbv^min = 100000, dataCN+1][6]; 

double qb=0,Al, B1, Cl, D1, Ta, Tb, Ttip, Tave, temp, qa, segment_diff; 

float error, V, epsilon, k, Tbase, density, mass, K2, eta, PN; 

int T, L, M, $, H, d, K, i, X, 2, Y; 

int overflow, possible_overflow = No, number_of_tries=0; 

void inputO, computeO, resultsO; 


jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


MAIN PROGRAM 


*****ifc**********i*******i**********************************i***********i*****:* k j 


int main(){ 

FILE *data_in, *out; 

data_in = fopenC'units.ext","r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&2,&Y,&error) 

fclose(data_in); 

error = error / 100; 

inputO; 

K2 = 2 * epsilon * sigma * v * pow(Tbase,4); 

Tave = Tbase * pow(V, 0,25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while (fabsCoverall_diff) > (Tbase * error) 11 i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * width * epsilon * pow(Ta, A) * depth; 

while (i <= N && overflow == No){ 

segment_diff = 10; 

while (fabs(segment_diff) > ((Tb * error) / N) && overflow == No){ 

A = A^O; B = B_(); C = C^O; D = D_(); 
temp =A*Ta+B* qa; 


77 


qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segmented!ff; 

if (Tb > Tbase + 50 j| Tb < 0)C 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv^min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;» 

if (overflow == No)C 

data[i][0] = qb; data[i][1] = Tb; data[i][2] = A; 
data[i][3] = B; data [i] [4] = C; data [i] [5] = D; 
Ta = Tb; qa = qb; ++i;>> 


if (number__of_tries > 29)<; 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(0);> 

if (overflow == No)<; 

if (Tmax < Ttip) 

if ( Tb < Tbase ) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overall_diff >0) Ttip=Ttip- inc; 
else Ttip = Ttip + inc;> 

else 

Ttip = Ttip - overall_diff; 

number_of_t ries++; 

if (Ttip < Tave) Ttip = Tave; 

gotoxy(27,24); 

cprintf("N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb));>> 
for (i = N - 1; i >= 1; 

A1 = A * data[i)[2] + B * data[i][4]; B1 = A * data[i][3] + B * data[i][5]; 

Cl = C * data[i][2] + D * data[i][4]; D1 = C * data[i) [3] + D * data[i][5]; 

A = A1; B = B1; C=C1; D=D1;> 

if (S == 1){ 

out = fopen(file_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf(out, "%e, %e, %e\n", w_L(d, i * length / N), w_T(T,data[i] [1]), w_H(H, 

data[i][0]) ); 

fclose(out);> 

Ta = Ttip; 

mass = width * length * depth * density * 1000; 

PN = length * sqrt( sigma * epsilon * pow(Tbase,3) / ( k * width)); 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

results(); 

exit(O); 

return;} 

y* *********************************************************************** 

OUTPUT PROCEDURE 

*★*★★*★★★★****★★*******★★★**★*****★★★**★***★★**★***★*****★****★**★★★****★* 


void results(void)C 


clrscrO; 

cprintf("\a THE LONGITUDINAL FIN OF RECTANGULAR PROFILE"); 

cprintf("\r\n ==================================================================\r\n"); 
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cprintfC*' Density = %A.2f kg / dm^3\r\n'\densi ty) ; 

cprintfC" Fin Depth = "); p_L(L, depth); 

cprintfC* Total Mass = '*); p_M(M, mass); 

cprintfC" Fin Length = "); p_L(L, length); 

cprintfC" Fin Width = "); p_LCd, width); 

cprintfC" Emissivity = %4.2f\r\n", epsilon); 

cprintfC" Profile Area = "); p_ACX, width*length); 

cprintfC" Profile Number = %4.3f\r\n",PN); 

cprintfC" Base Heat Flow = "); p_HCH, qb); 

cprintfC" Fin Efficiency = %4.3f\r\n",eta); 

cprintfC" Base Temperature = "); p_TCT,Tb); 

cprintfC" Edge Temperature = "); p_TCT,Ta); 

cprintfC" Thermal Conductivity = "); p_KCK, k); 

cprintfC" Environmental Parameter = %4.2f\r\n", V); 

cprintfC" ================rr==r===================r===========rr============\r\n"); 

cprintfC" Transmission Parameter Matrix \r\n\r\n"); 

cprintfC" A = %4.4f B = ", A); p_ZCZ,B); 

cprintfC" C = "); p_YCY,C); 

cprintfC" D = %4.4f", D); 

getchc); 

return;} 


j-k *********************************************************************** 

INPUT PROCEDURE 

kkkkkkkkkkkkkkkkkkkkkk-k'kkk'k'kkk-kkkkkit-k-k-kk-k'kk'k'k-k-kk-k-k-kkkk'kkkkkk'k'k'kk'k'iftfk-k-kkkk-k k j 


void inputCvoidX; 

cIrscrC);_setcursortypeC_NORMALCURSOR); 


cprintfC"\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROFILE\n\r"); 

cprintfC"\r\n ==============================^================r==================\r\n"); 

cprintfC" Fin Depth "); depth = i_LCL); 

cprintfC" Fin Length "); length = i_LCL); 

cprintfC" Fin Width "); width = i_LCd); 

cprintfC" Emissivity = "); scanfC"%f", Sepsilon); 

cprintfC" Density C kg / dm^3 ) = "); scanfC"%f", Sdensity); 

cprintfC" Base Temperature "); Tbase = i_TCT); 

cprintfC" Thermal Conductivity "); k = i_KCK); 

cprintfC" Environmental Parameter CO to .8) = "); scanfC"%f", &V); 

if CS == 1) C 

cprintfC" Output Filename is = "); 


scanfC"%s", file_name);> 


cprintfC"\r\n=================================================================="); 

if C epsilon > 1.0 || epsilon <= 0)<; 

cprintfC"\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchC); exitCO);} 

if C V < 0 I 1 V > 0.8)<: 

cprintfC"\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 
getchc); exitCO);} 

_setcursortypeC_NOCURSOR); 
return;} 

jk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

GLOBAL FUNCTIONS 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk k j 


double hrCvoid)C 

float Tav = CTb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) - CK2 j * Tav)));} 
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double m(void)<; 

return( sqrt( ( 2.0 * hr()) / ( k * width) ) );> 


double Yb(void){ 

return( sqrt( 2.0 * hr() * k 

double A_(void)<; 

returnC cosh( m() length / 

double B_(void){ 

return ( sinh( m() * length / 

double C_(void){ 

return ( sinh( m() * length / 

double D_(void)<; 

return( cosh( m() length / 


* width) * depth );> 
N ) );> 

N ) / Yb() );> 

N ) * YbO );> 

N ) );> 
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j-k ★***★★★★★★★*********★***★*★★★*******************★*★********★★★★*******★**** 


PERFORMANCE Analysis Program 
Longitudinal Fin of Trapezoidal Profile 

*****4r**************************4r*******************************i«r*********** k j 


^include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <conio.h> 
#include "bessel.h'* 
#include "units.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 100 

sigma 5.66961e-8 


char file_name[20]; 

double A_(), B_(), C_(), D_(), m(), hr(), uCdouble x), A, B, C, D; 

double length, depth, edge_width, base_width, fin_taper; 

double b1, db, a, b, ua, ub, Ta, Tb, a_width, b_width; 

double qb=0, inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data[N+1][6] 

double A1, B1, Cl, Dl, taper_ratio, Ttip, temp, qa, segment_diff; 

float V, epsilon, k, error, Tbase, density, mass, K2, Ap, eta, PN; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

int overflow, possible_overflow = No, number_of_tries=0; 

void inputO, resultsO; 


^k kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


MAIN PROGRAM 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk kj 


int main(){ 

FILE *data_in, *out; 

data_in = fopenC'units.ext","r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error = error / 100; 

inputO; 

taper_ratio = edge^width / base_width; 
b1 = length / ( 1 ” taper_ratio ); 
fin^taper = atan2( base_width, 2*b1 ); 
db = length / N; 

K2 = 2 * epsilon * sigma * v * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrtC 2 * V )))/2); 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while (fabsCoveral l_diff) > (Tbase * error) 1 j i <= N )<; 


overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * length; 

a = b1 - length; 
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while (i <= N && overflow == No)( 


segment_dfff = 10; 

a_width = edge^width a / ( b1 - length); 
b_width = edge_width ’^(a + db)/(b1- length); 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No)( 

ua = u(a) ; ub = u(a + db); 

A = A_(); B = B_(); C = C_(); D = D_(); 
temp = A*Ta + B’^qa; 
qb = C * Ta + D * qa; 
segmented!ff = ( temp - Tb) / 2; 

. Tb = temp + segment_diff; 

if (Tb > Tbase +50 || Tb < 0){ 

overflow = Yes; possible_overfIow = Yes; 
if (Ttip < Tbv^min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 


if (overflow == No){ 

datati][0] = qb; data[i][1] = Tb; data[i] [2] = A; 
data[i][3] = B; datati][4] = C; data[i] [5] = D; 
Ta = Tb; qa = qb; a = a + db; ++i;}} 


if (number_of_tries > 29){ 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(O);} 

if (overflow == No){ 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb ^ Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;> 

else 

Ttip = Ttip - overall_diff; 
number_of_t ries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i-OC 

A1 = A ^ datati]t2] + B data[i]t4]; Bl = A * data[i][3] + B data[i][5]; 

Cl = C datati]t2] + D data[i]t4]; D1 = C data[i][3] + D * data[i]t5]; 

A = A1; B = Bl; C=C1; D=D1;} 

if (S == 1){ 

out = fopen(file_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf(out, "%e, %e, %e\n", w_L(d, i * length / N), w_T(T,data[i] tl]), w_H(H, 

dataCi] [0]) ); 

fclose(out);} 


Ta = Ttip; 

Ap = .5 length ( edge_width + base_width ); 
mass = Ap * length * density * 1000; 

PN = 2 sigma * epsilon * pow( length, 2) * pow(Tbase,3) / (k* base_width); 
eta = qb / ( 2 sigma * epsilon * length * length * pow(Tbase,4)); 
results(); 
exi t(0); 
return;} 
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OUTPUT PROCEDURE 

★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★********<k*******<k<t*********<k******** -k j 


void results(void){ 


clrscrC); 
cprintf(" 
cprintf("\r\n 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
getch(); 
return;> 


THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 


Density 
Fin Depth 
Total Mass 
Fin Length 
Emissivity 
Profile Area 
Profile Number 
Base Heat Flow 
Fin Efficiency 
Fin Edge Width 
Fin Base Width 
Base Temperature 
Edge Temperature 
Thermal Conductivity 
Environmental Parameter 




%4.2f kg / dm^3\r\n",density); 
"); p_L(L, depth); 

"); p M(M, mass); 

"); pIlcL, length); 

%4.2f\r\n", epsilon); 

"); p_A(X,Ap); 

%4.3f\r\n",PN); 

"); p_H(H, qb); 

%4.3f\r\n",eta); 

"); p_L(d, edge_width); 

"); p_L(d, base_width); 

"); p_T(TJb); 

"); p T(T,Ta); 

"); p_K(K, k); 

%4.2f\r\n", V); 


Transmission Parameter Matrix \r\n\r\n"); 
A = %4.4f B = ", A); p 2(2,B); 

C = "); p_Y(Y,C); 

D = %4.4f", D); 


=\r\n"); 






INPUT PROCEDURE 


kkkkkkkitk'k 


kkkkkk-kkk-kk-kk 


kkkkkkkkkkk-kkkkk 


kkkkkkkifkkik 


void input(void){ 

clrscrC); _setcursortype(__NORMALCURSOR); 


cprintf("\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE\r\n"); 

cpr intf ("\r\n =======r===r===r==r===rrr=====r============================r======\r\n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Fin Length "); length = i_L(L); 

cprintfC Fin Edge Width "); edge_width = i_L(d); 

cprintfC Fin Base Width "); base_width = i_L(d); 

cprintfC Emissivity = "); scanfC%f", &epsilon); 

cprintfC Density ( kg / dm^3 ) = "); scanfC%f", &density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanfC%f", &V); 

if (S == 1) f 

cprintfC Output Filename is = "); 


scanfC%s", file_name);> 


cprintfC\r\n =================================r==r=================r==========="); 

if ( edge_width >= base__width)( 

cprintfC\r\n\r\nFATAL ERROR: Fin Edge Width >= Fin Base Width"); 
getchO; exit(0);> 

if ( epsilon > 1.0 [j epsilon <= 0){ 

cprintfC\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchO; exit(0);> 
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if ( V < 0 I I V > 0.8)C 

cprintf("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8”) 
getchO; exit(O);} 

_setcursortype(_NOCURSOR); 
return;} 

j-k kkkkkkkilickkitkk-kk-kkkkkkkkkkkkkkkkkkkkkkkkifkkifkk-kk-kkkkkkkkkkkkkkkkkifkkkkk 

GLOBAL FUNCTIONS 

kkkkkkkkkkkkkkkkkkkkkkkkk************************************************* * j 


double hr(void)<; 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - (K2 / (2 Tav)));> 

double m(void)<; 

returnC hr() / (k * sin(fin_taper)));> 

double uCdouble x)<; 

returnC 2 * sqrtCx * mC) ) );> 

double A_Cvoid){ 

returnC ua * C11Cua)*K0Cub) + l0Cub)*K1Cua)) );> 

double B_Cvoid)<; 

return C ub’*'ub’*'CI0Cub)*K0Cua)- I0Cua)*K0Cub)) / C 2*k*b_width*depth*mC)));> 
double C_Cvoid){ 

return C 2’*k’*a_width*depth’*'mC)’*'ub*C 11 Cub)*K1 Cua) - 11 Cua)’^K1 Cub))/Cua) );> 
double D_Cvoid)<; 

returnC ub*C l0Cua)*K1Cub)+I1Cub)*K0Cua)) );> 
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y* ★★★**★**★***★**★★**★**★★*★**★★★★★★*****★*★★**★**★★★★★★★★★*★**★★★★**★***★★★★ 


PERFORMANCE Analysis Program 
Longitudinal Fin of Triangular Profile 


#include <stdio.h> 
#include <stdlib.h> 
^include <math.h> 

#include <conio.h> 
#include "bessel.h" 
#include "units.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 99 

sigma 5.66961 e-8 


char file_namet20]; 

double A_(), B_(), C_(), D_(), m(), hr(), uCdouble x), A, B, C, D; 

double length, depth, edge_width, base_width, fin_taper; 

double b1, db, a, ua, ub, Ta, Tb, a_width, b_width; 

double qb=0, inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, datatN+13[6] 

double A1, B1, Cl, D1, Ttip, temp, qa, segment_diff; 

float V, epsilon, k, error, Tbase, density, mass, K2, Ap, PN, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

int overflow, possible_overflow = No, number_of_tries=0; 

void inputO, resultsO; 


y* **★★*****★★**★★****★*********★★★***★★★★★★*****★★****★*★*★****★★******★★ 


MAIN PROGRAM 


*******'*itf^'*'*:t^****i*i****'**********i*i*******'*************4r***ir4r**'****4r4rir*iAr'4ririririr ★ j 


int mainOC 

FILE *data_in, ^out; 

data_in = fopenC'uni ts.ext“,'‘r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f•',&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&er^o^); 

fclose(data_in); 

er^or = error / 100; 

inputC); 

db = length / ( N + 1 ); 

b1 = length; 

length = length - db; 

fin_taper = atan2( base_width, 2*b1 ); 

K2 = 2 * epsilon * sigma ^ V ^ pow(Tbase,4); 

Tave = Tbase * powCV, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while (fabsCoverall_diff) > ( Tbase * error) 1[ i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segment_diff = 10; 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No){ 
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ub = u(db); 

temp = Ta / lO(ub); 

C = 2 * k * edge_width * depth * pow(m(),2) * lUub) / (ub * lO(ub)); 
qb = C * temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase +50 ]| Tb < 0){ 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min * inc); 
if (Ttip < lave) Ttip = Tave;>> 

if (overflow == No)<! 

data[i] [0] = qb; data[i][l] = Tb; data[i][2] = 1; 
data[i][3] = 0; data[i][4] = C; data[i][5] = 1; 

Ta = Tb; qa = qb; ++i;> 


a = bl ‘ length; 

while (i <= N && overflow == No){ 
segment_diff = 10; 

a_width = edge_width * a / ( bl - length); 
b_width = edge_width * ( a + db ) / ( bl - length); 

while (fabs(segment_di ff) > ( Tb * error) / N && overflow == No)<; 

ua = u(a) ; ub = u(a + db); 

A = A_(); B = B_(); C = C_(); D = D_(); 
temp =A*Ta+B* qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase +50 11 Tb < 0){ 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 

if (overflow == No){ 

data[i][0] = qb; data[i][l] = Tb; data[i][2] = A; 
data[i][3] = B; data[i][4] = C; data[i][5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;>> 


if (number_of_tries > 29){ 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(0);> 

if (overflow == No)<; 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overall_diff >0) Ttip = Ttip* inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overall_diff; 
number_of_tr i es++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 
if (Ttip < Tave) Ttip = Tave;}} 
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for (i = N - 1; i >= 1; i--)C 

A1 = A * data[i][2] + B data[i3[4]; B1 = A * data[i][3] + B * data[i][53; 

Cl = C * dataCi3C23 + D * data[i3[43; D1 = C * data[i3 [33 + D * data[i3[53; 

A = Al; B = Bl; C=C1; D=D1;> 

length = length + db; 

if (S == DC 

out = fopenCfile_name, "w”); 
for (i = 1; i <= N; i++) 

fprintfCout, '‘%e, %e, %e\n", w_L(d, i * length / N), w_T(T,data[i3 [13), w H(H, 

data[i3[03)); 

fcloseCout);> 

Ta = Ttip; 

Ap = .5 * length * base_width; 
mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow(length,2) * pow(Tbase,3) / (k*base_width); 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

results(); 

exit(O); 

return;} 

OUTPUT PROCEDURE 

★ A************************************************************************ Tk y 


void results(void)C 


cl rscrC); 

cprintf("\a THE LONGITUDINAL FIN OF TRIANGULAR PROF ILE\r\n"); 

cprintf(" =====r=:===========================================================\r\n“); 

cprintfC Density = %4.2f kg / dm^3\r\n" ,densi ty); 

cprintfC Fin Depth = *'); p_L(L, depth); 

cprintf(“ Total Mass = "); p_M(M, mass); 

cprintf(“ Fin Length = "); p_L(L, length); 

cprintfC Emissivity = %4.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X,Ap); 

cprintfC Profile Number = %4.3f\r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4.3f\r\n“,eta); 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC Base Temperature = "); p_T(T,Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC ==================================================================\r\n"); 

cprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = %4.4f B = ", A); p^Z(Z,B); 

cprintfC C = "); p_Y(Y,C); 

cprintfC D = %4.4f'‘, D); 

getchO; 

return;} 


Ilr'****'*'***********'***********4r4r4r*4r4r4r4r4r4r4r*4r4r4r4n*n*r4r4r4r4r4r4r4r4r4r4r4r4r4r4r******4r4r*** 

INPUT PROCEDURE 

*';Ar'4r'4r4r4r*ilr4r4r'*'*ir**'A'*'*'A'4r'<r4r4r*4r4r*'*'****4r4r4r4r4r'4r4r*4r4r4r4r*'*'4r'*'*ifr******4r4r'*'*'4r4r*4r4r4r'4r'*'4rir4r4r'4r4r ★ j 


void input(voidx; 

clrscrC); _setcursortype(_NORMALCURSOR); 


cprintfC\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROF I LE\r\n"); 

cprintf("\r\n ==================================================================\r\n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Fin Length "); length = i_L(L); 
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cprintfC" Fin Base Width "); base_width = i_L(d); 

cprintfC" Emissivity = "); scanf("%f", fiepsilon); 

cprintfC" Density ( kg / dm^3 ) = "); scant^density); 

cprintfC" Base Temperature "); Tbase = i_TCT); 

cprintfC" Thermal Conductivity "); k = i_KCK); 

cprintfC" Environmental Parameter CO to .8) = "); scanfC"%f", &V); 

if CS == 1) ( 

cprintfC* Output Filename is = "); 

scantC"%s", file_name);} 

CprintfC"\r\n = = = = = = = = = = = = r = === = = = rrr = = = r = = ===r = = = = = = = = = = = = = = === = === = r = = = = = = = = r = «'); 


if C epsilon > 1.0 jj epsilon <= 0)( 

cprintfC"\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchC); exitCO);} 

if C V < 0 11 V > 0.8){ 

cprintfC"\r\n\r\nFATAL ERROR; Environmental Parameter is not between 0 and 0.8"); 
getchC); exit CO);} 

edge_width = base_width / le2; 

_setcursortypeC_NOCURSOR); 
return;} 

j-k 'k'k'k'k'kkk-k-kkk'k'k-k'kifkiek-k-kieickkk-k-k-k'k-kifkifickk-tc-k-k-kickickickk-k'k'k'k-k-k'k-kie-kkkick-tck'k-k-k'k-k-k'k 

GLOBAL FUNCTIONS 

■kk'k-k-k-kkkkkkkkk-kkkkkkkkk-kkkkk-k-kkkkit kkkk'k'kk-k-k kkk'kkk-k kk-kk-kkifkifk kick •kk'k'k'kifk kick k k j 


double hrCvoid)( 

float Tav = CTb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) - CK2 / C2 * Tav)));} 
double mCvoid)( 

returnC hrC) / Ck * sinCfin_taper)));} 

double uCdouble x)<; 

returnC 2 * sqrtCx * mC) ) );} 

double A_Cvoid)( 

returnC ua * C11Cua)*K0Cub) + I0Cub)*KlCua)) );} 

double B_Cvoid)( 

return C ub*ub*CI0Cub)*K0Cua)-I0Cua)*K0Cub)) / C 2*k*b_width*depth*mC)));} 
double C_Cvoid)( 

return C 2*k*a_width*depth*mC)*ub*C11Cub)*KlCua) - 11Cua)*KlCub))/Cua) );} 
double D_Cvoid)C 

returnC ub*C I0Cua)*KlCub) + nCub)*K0Cua)) );} 
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/ 




A******* 


★ A* AAAAAAAAAAA 


AAAAAAAAAAAAAAA 


PROFILE AREA Analysis Program 
Longitudinal Fin of Rectangular Profile 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa a j 


^include <stdio.h> 
^include <stdlib.h> 
#include <math.h> 

#include <conio.h> 

#include "units.h" 


#define Yes 1 

#define No 0 

^define N 100 

#define sigma 5.66961e-8 


char file_name[20] ; 

double length, width, depth, starting_width, ending_width, width_inc; 

double A, B, C, D, Ta, Tb, qb=0; 

float V, IC2, epsilon, k, Tbase, density, mass, Ap, error, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

void inputO, computeO; 


J'k aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


MAIN PROGRAM 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA -k j 


int mainOC 

FILE *data_in, *out; 

data_in = fopenC'units.ext","r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 
fclose(data_in); 
error = error / 100; 

inputO; 

width^inc = (ending_width - starting_width ) / 100; 
width = starting_width; 

while (width < ending_width)<; 

compute(); 

while (qb == OX 

width = 1.001 * width; 
computeC);} 

gotoxyC19,23); 

cprintf("width = %4.1f length = %A.1f heat = %4.1f", w_L(d,width), w L(L,length), 
w_H(H, qb)); 

if (S == DC 

out = fopen(file_name, "a"); 

fprintf(out, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n", w_L(d,width), 

w_T(T,Tb), w_T(T,Ta), eta, w_H(H, qb), w_M(M,mass), A, w_,Z(Z,B), w_Y(Y,C), D) 
fclose(out);) 

width = width + width^inc;} 


gotoxy(21,2A);clreoK);cprintf("\aCompleted \r\n"); 
clreol(); 
getch(); 
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exit(O); 
return 0;} 

jit •k’k'k-k'k-k-k-k-k-tt-k-k-kit'k-k-k-k-k-k-k'k'k-k'k-k'k-k-k-k-k-k-k-k-itit-k'kickifk-k-k'tfk'ifkickick'k-k'tfk-k'k-ictftctck-k-kickifk'k-k 

COMPUTE ROUTINE 

****★*****★★*★**★★*★★★**★★**★*★★*★*★★***★****★**★*★★★★★*★★★★★**★*★★*★*★*** -k j 


void compute(void){ 

double A_(), B_(), CJ), D_(), m(), hr(), Yb(); 

double inc, overall^diff = 100, Tmax = 0, Tbv_min = 100000, data[N+1][6]; 

double Al, B1, Cl, D1, Ttip, temp, qa, Tave, segmented]ff; 

int overflow, possible_overflow = No, number_of_tries=0; 

length = Ap / width; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2) 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while (fabsCoveral l_di f f) > ( Tbase * error) | 1 i <= N )<; 

overflow = No; i =1; Ta = Ttip; Tb = Ta; 

qa = sigma * width * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == No){ 

segment_diff = 10; 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No)< 

A = A_(); B = B_(); C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase + 50 11 Tb < 0){ 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No)C 

data[i][0] = qb; data[i][1] = Tb; data[i)[2] = A; 
data[i]C3) = B; data[i3[4] = C; data(i] [5] = 0; 

Ta = Tb; qa = qb; ++i;}} 

if (number_of_tries > 29){ 
qb=0;return;} 

if (overflow == No){ 

if (Tmax < Ttip) 

if (Tb < Tbase ) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overall^diff >0) Ttip=Ttip' inc; 
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else 


else Ttip = It ip + inc;} 


Ttip = Ttip “ overall_diff; 
number_of_t ries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 

if (Ttip < Tave) Ttip = Tave;}} 

for (i = N " 1; i >= 1; i “ )<; 

A1 = A * data[i][2] + B * data[i)[4]; B1 = A * data[i] [3] + B * data[i)[5]; 

Cl = C * data[i) [2] + D * data[i) [4]; D1 = C * data[i] [3] + D * data[i) C5] ; 

A = Al; B = Bl; C=C1; D=D1;} 

Ta = Ttip; 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

mass = width * length * depth * density * 1000; 

return;} 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkifkkkkkkkit'kkifkifkifkkkkkkkitkitkkkkkitititkickkitkkit 

INPUT PROCEDURE 


void input(voidx; 

clrscrC); _setcursortype(_NORMALCURSOR); 


cprintf("\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROFILE"); 

cprintf("\r\n ==================================================================\r\n"); 

cprintfC" Fin Depth "); depth = i_L(L); 

cprintfC" Starting Fin Width "); starting_width = i_L(d); 

cprintfC" Ending Fin Width "); ending_width = i_L(d); 

cprintfC" Profile Area "); Ap = i_A(X); 

cprintfC" Emissivity = "); scanfC"%f", Sepsilon); 

cprintfC" Density C kg / dm^3 ) = "); scanf("%f", Sdensity); 

cprintfC" Base Temperature "); Tbase = i_T(T); 

cprintfC" Thermal Conductivity "); k = i_lCCK); 

cprintfC" Environmental Parameter (0 to .8) = "); scanfC"%f", &V); 

if (S == 1) C 

cprintfC" Output Filename is = "); scanf("%s", file_name);} 

cprintfC"\n ==================================================================\r\n"); 


if C epsiI on > 1.0 j | epsiI on <= 0)C 

cprintf("\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchC); exit(O);} 

if C V < 0 1 1 V > 0.8)<: 

cprintf("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 
getchC); exitCO);} 

_setcursortype(_NOCURSOR); 
return;} 

y* kkifickk-kk* ******************** kk-kkkk-kitkkkkkkkkkkifickk-kkkicifififkkit ********** 

GLOBAL FUNCTIONS 

************************************************************************** *! 


double hrCvoidX; 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * powCTav, 3) • C)C2 / (2 Tav)));} 

double mCvoidX 

returnC sqrtC C 2.0 * hrO) / C k * width) ) );} 
double YbCvoidX; 
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returnC sqrt( 2.0 * hr() * k * width) * depth );> 


double A_(void){ 

returnC cosh( m() * length / N 

double B_(void)< 

return ( sinhC m() length / N 

double C_(void)<: 

return ( sinh( m() * length / N 

double D_(void)<; 

returnC coshC mC) * length / N 


);> 

) / YbC) );> 
) * YbC) );> 
);> 
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I* *★★*****★★★*****★★*****★*★***★**********★*★*★★*** 




PROFILE AREA Analysis Program 
Longitudinal Fin of Trapezoidal Profile 

ititit'kitit'k'k-k-kif-k-k-k-k'k-k'k'k-k'k'k'k-k-k-k-k-k-k-kit-k-k-k-k-k-it-tt'k'kitit-k’kit-k-k-k-k-k'k'k'kit-kitit'kit'k'k-k-kit-k'k'k-k'k'kitit-kit-kit ★ j 


#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 

#include <conio.h> 
^include “units.h*' 
#include “bessel.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 100 

sigma 5.66961e-8 


char f i le_name[20]; 

double length, edge_width, base_width, a_width, b__width, depth, starting_width, ending_width 

double A, B, C, D, Ta, Tb, qb=0, ua, ub, fin_taper, a, width_inc; 

float V, K2, epsilon, k, Tbase, taper_ratio, density, mass, Ap, error, eta; 

int I, L, M, S, H, d, K, i, X, 2, Y; 

void inputO, computet); 




MAIN PROGRAM 




int main(){ 

FILE *data__in, *out; 

data_in = fopenC'uni ts.ext'',"r''); 

fscanf(data_in,“%d %d %d %d %d %d %d %d %d %d %f'',&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 

fclose(data_in); 

error = error / 100; 

input(); 

width_inc = (ending_width - starting_width ) / 100; 
base_width = starting_width; 

while (base_width < ending_width)( 

computet); 

while tqb == 0){ 

base_width = 1.001 ^ base_width; 
computet);> 

gotoxyt19,23); 

cprintff'width = %4.1f length = %4.1f heat = %4.1f'', w_Ltd,base_width), w_LtL, length), 
w_HtH, qb)); 

if tS == 1){ 

out = fopentfile_name, “a“); 

fprintftout, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n“, w Ltd,base_width), 

w_TtT,Tb), w_TtT,Ta), eta, w^HtH, qb), w__MtM,mass), A, w__Zt2,'B), w__YtY,C), D); 
fclosetout);> 

base_width = base_width + width__inc;> 


gotoxyt21,24);cl reol t) ;cprintft"\aCompleted \r\n''); 
clreolt); 
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getchO; 
exit(0); 
return 0;> 

jit **★****★***********★**★****★★★★*********★★★★****★**★★★★★★**★★*★****★★★★ 

COMPUTE ROUTINE 

************************************************************************** 


void cofnpute(void){ 

double A_(), B_(), C_(), D_(), m(), hr(), uCdouble x); 

double length_inc, b1, db; 

double inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, data[N+1][6]; 

double A1, B1, Cl, D1, Ttip, temp, qa, segment_diff; 

fnt overflow, possible_overflow = No, number_of_tries=0; 

edge_width = base_width * taper_ratio; 
length = 2 * Ap / ( base_width + edge_width ); 
b1 = length / ( 1 - taper_ratio ); 
fin_taper = atan2( base_width, 2*b1 ); 
db = length / N; 

K2 = 2 * epsilon * sigma * v * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave 0.0001; 
inc = .5 * Ttip; i=1; 

while (fabsCoverall_diff) > ( Tbase * error) 1] i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = b1 - length; 

while (i <= N && overflow == No){ 
segment_diff = 10; 

a_width = edge_width * a / ( b1 - length); 
b_width = edge_width * ( a + db ) / ( b1 - length); 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No)( 

ua = u(a) ; ub = u(a + db); 

A = A_(); B = B_(); C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb =C*Ta+D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase + 50 ] ] Tb < 0)<: 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv^min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 

if (overflow == No)<; 

data[i][0] = qb; data[i][13 = Tb; data[i3[2] = A; 
data[i][3] = B; data[i][4] = C; data[i3[53 = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}> 


if (number_of_tries > 29)C 
qb=0; return;) 
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if (overflow == No)( 

if (Tmax < Trip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes)( 
inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overall_diff; 
number_of _t r i es-*’-*'; 
gotoxy'(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 

if (Ttip < Tave) Ttip = Tave;>> 

for (i = N - 1; i >= 1; i--)( 

A1 = A * data[i][2] + B * data[i][4]; B1 = A * data[i][3] + B * data[i][5]; 

Cl = C * data[i][2] D * data[i][4]; D1 = C * datali] [3] + D * data[i][5]; 

A = A1; B = B1; C=C1; D=Dl;} 

Ta = Ttip; 

mass = .5 * ( edge_width + base^width ) * length * depth * density * 1000; 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase, 4)); 
return;} 

INPUT PROCEDURE 

★ ★*★★*★★★*★★★★★★★******★***★★*★*★★★★*****★★★★★★★★*★★★**★★**★★★**★★★**★★*★* it j 


void input(void)( 

clrscr(); _setcursortype(_NORMALCURSOR); 


cprintfC THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 

cprintf("\r\n ==================================================================\r\n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Starting Base Fin Width "); starting_width = i_L(d); 

cprintfC Ending Base Fin Width "); ending_width = i_L(d); 

cprintfC Taper Ratio = "); scanfC%f‘', &taper_rat io); 

cprintfC Profile Area "); Ap = i_A(X); 

cprintfC Emissivity = "); scanfC%f", &epsilon); 

cprintfC Density ( kg / dm^3 ) = "); scanfC%f", &density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanfC%f", &V); 

if (S == 1) ( 

cprintfC Output Filename is = "); scanfC%s", file_name);} 

cprintfC\n ===========================================================r======\r\n"); 


if ( taper_ratio > .99 |] taper_ratio <= 0.0) ( 

cprintfC\r\n\r\nFATAL ERROR: Taper Ratio is not between 0 and 1"); 
getchO; exit(O);} 

if ( starting^width >= ending_width)( 

cprintf("\r\n\r\nFATAL ERROR: Starting Base Fin Width >= Ending Base Fin Width"); 
getchO; exit(O);} 

if ( epsilon > 1.0 || epsilon <= 0){ 

cprintfC\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchO; exit(O);} 

if ( V < 0 1[ V > 0.8){ 

cprintfC\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 
getchO; exit(O);} 
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_setcursortype(_NOCUR$OR); 
return;) 

l^t A********************************************************************** 

GLOBAL FUNCTIONS 

************************************************************************** *! 


double hr(void){ 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));> 
double m(void)<; 

returnC hr() / (k * sin(fin_taper)));} 

double uCdouble x)<; 

returnC 2 * sqrtCx * m() ) );> 

double A_(void){ 

returnC ua * C11Cua)*K0Cub) + I0Cub)*KlCua)) );} 

double B_Cvoid){ 

return C ub*ub*ClOCub)*KOCua)-lOCua)*KOCub)) / C 2*k*b_width*depth*mC)));} 
double C_Cvoid){ 

return C 2*k*a_width*depth*mC)*ub*C11Cub)*K1Cua) • 11Cua)*K1Cub))/Cua) );} 
double D_Cvoid){ 

returnC ub*C 10Cua)*K1Cub)+11Cub)*K0Cua)) );} 
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j* *************************************************************************** 


PROFILE AREA Analysis Program 
Longitudinal Fin of Triangular Profile 

'kii-kit-k-k-k-k-k-k-k-k-kit-k-k-k-k-k-k-k-k-k-kii-k-k-k-kit-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k ★ j 


#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <conio.h> 
#include "units.h" 
#include "bessel.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 99 

sigma 5.66961e-8 


char file_name [20]; 

double length, width_edge, base_width, a_width, b_width, depth, starting^width, ending_width 

double A, B, C, D, Ta, Tb, qb=0, ua, ub, fin_taper, a, inc_width; 

float V, K2, epsilon, k, Tbase, density, mass, Ap, error, eta; 

int T, L, M, S, H, d, K, i, X, Z, Y; 

void inputO, computeC); 


j-k ******★★★**★★***★★***★★★★★******■*★*'***★*★**★★★**★★**★★★★*★★*★*★★★★★★★★★ 


MAIN PROGRAM 


**★★★***★**★**★**★★******★★★**★****★★★★***★★*★★★***★**★*★★★**★★★★*★★★*★★★* -k j 


int rnainCX 

FILE *data_in, *out; 

data_in = fopenC'units.ext","r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&Z,&Y,&error); 
fclose(data_in); 
error = error / 100; 

input(); 

inc_width = (ending_width - starting_width ) / 100; 
base_width = starting_width; 

while (base_width < ending_width){ 

computeC); 

while (qb == 0){ 

base_width = 1.001 * base_width; 
computeC);) 

gotoxyC19,23); cprintfC'width = %4.1f length = %4.1f heat = %4.1f", w_LCd,base_width), 
w_LCL,length), w_HCH, qb)); 

if CS == 1){ 

out = fopenCfile_name, "a"); 

fprintfCout, "%e, %e, %e, %e, %e, %e, %e, %e, %e, %e\n", w_LCd,base_width), 

w_TCT,Tb), w_TCT,Ta), eta, w_HCH, qb), w_MCM,mass), A, w_ZCZ,B), w_YCY,C), D); 
fcloseCout);> 

base_width = base_width + inc_width;> 


gotoxyC21,24);clreolC);cprintfC"\aCompleted \r\n"); 

clreolC); 

getchC); 
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exit(0); 
return 0;> 

j-k ★*★*★***★★★★***★★★★***********★★**★*******★*★******★***★*****★★*★★*★★** 

COMPUTE ROUTINE 


void compute(void){ 

double A_(), B_(), C^O, D^O. m(), hr(), u(double x); 

double length_inc, b1, db; 

double inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000, dataCN+l][6] 

double Al, Bl, Cl, Dl, Ttip, temp, qb, segment_diff; 

int overflow, possible_overflow = No, number_of_tries=0; 

length = 2 * Ap / base_width; 

width_edge = base_width / 100; 

db = length / ( N + 1 ); 

bl = length; 

length = length ■ db; 

fin_taper = atan2( base_width, 2*b1 ); 

K2 = 2 * epsilon * sigma * V * pow(Tbase,A); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2) 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while (fabsCoveralL_diff) > ( Tbase * error) ]] i <= N ){ 

overflow = No; i =1; Ta = Ttip; Tb=Ta; 
segment_diff = 10; 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No){ 

ub = u(db); 

temp = Ta / lO(ub); 

C = 2 * k * width_edge * depth pow(n(),2) * Il(ub) / (ub * lO(ub)) 
qb = C * Tb; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase +50 | | Tb < 0)<; 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 

if (overflow == No)<; 

data[i]C0] = qb; data[i][13 = Tb; data[i] [2] = 1 ; 
data[i][3] = 0; data[i][4] = C; dataCi]C5] = 1; 
a = bl - length; 

Ta = Tb; qa = qb; ++i;> 

while (i <= N && overflow == No){ 

segment_diff = 10; 

a_width = width_edge * a / ( bl • length); 
b^width = width_edge * ( a + db ) / ( bl - length); 

while (fabs(segment_di f f) > ( Tb * error) / N && overflow == No)(; 
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ua = u(a) ; ub = u(a + db); 

A = A_(); B = C = C_(); D = D_(); 

temp =A*Ta+B* qa; 
qb =C*Ta+D* qa; 
segmented!ff = ( temp - Tb) / 2; 

Tb = temp segment_di f f ; 

if (Tb > Tbase + 50 | | Tb < 0)<: 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv^min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 


if (overflow == No)( 

data[i][0] = qb; data[i][1) = Tb; data[i][2] = A; 
data[i][3) = B; data[i][4] = C; data[i][5] = D; 
Ta = Tb; qa = qb; a = a + db; ++i;>> 

if (number_of_tr 1 es > 29){ 
qb=0; return;} 

if (overflow == No)< 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 


overall_diff = (Tb - Tbase) / 2; 

if ( possible_overf low == Yes)<; 

inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;> 

else 

Ttip = Ttip - overall_diff; 
number_of_t ries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 

if (Ttip < Tave) Ttip = Tave;>> 


for (i = N • 1; i >= 1; i-OC 

AT = A * dataCi] [2] + B * dataCi] [4]; B1 = A * dataCi] [3] + B * dataCi] [5] ; 
Cl = C * dataCi) [2] + D * dataCi] [4]; D1 = C * dataCi) C3) + D * dataCi)C5); 
A = AT; B = Bl; C=C1; D=D1;} 


length = length + db; 

Ta = Ttip; 

mass = .5 * base_width * length * depth * density * 1000; 

eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

return;} 


j-k ★★★**★**★★♦*****★***★★★★*■*★★★★**★★★★★★★★**★★★★★★★**★*****★**★★**★★★*★** 


INPUT PROCEDURE 

'k-k-k-k-kk-k-k-k-k-k-k-k-kkkkk-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'kkk-kk-kkkk-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-kk-k-k-k-kictfk-k * j 


void input(void)<; 

cIrscr();_setcursortype(_NORMALCURSOR); 


cprintfC THE LONGITUDINAL FIN OF TRIANGULAR PROFILE\r\n"); 

cprintf(" =========================r========================================\r\n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Starting Base Fin Width *'); starting_width = i_L(d); 

cprintfC Ending Base Fin Width "); ending_width = i_L(d); 

cprintfC* Profile Area *'); Ap = i_A(X); 

cprintfC Emissivity = "); scanf("%f", &epsilon); 

cprintfC* Density ( kg / dm^3 ) = *'); scanf("%f", &density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC* Thermal Conductivity “); k = i_K(K); 
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cprintf( 


Environmental Parameter (0 to .8) = "); scanf("%f", &V); 


cprintfC Output Filename is = "); scanf("%s", file_name);) 

cprintf("\n ====================r=================rr==rr=rs=ssssss=s==r=======\r\n"); 


if ( starting_width >= ending^widthX 

cprintf("\r\n\r\nFATAL ERROR: Starting Base Fin Width >= Ending Base Fin Width") 
getchO; exit(O);} 

if ( epsilon > 1.0 jj epsilon <= 0)<; 

cprintf("\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchO; exit(0);> 

if ( V < 0 11 V > 0.8)C 

cprintf("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8") 
getchO; exit(0);> 

_setcursortype(_NOCURSOR) ; 
return;} 

*i(if-kit-k-k'ki(**-kif*-k*-k*****-k-k-k-k****ifk-k*-kic*(***i(***-k**************ififk***-k-k**** 

GLOBAL FUNCTIONS 


double hr(void){ 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));} 
double m(void)C 

return( hr() / (k * sin(fin_taper)));} 

double u(double x)<; 

return( 2 * sqrt(x * m() ) );} 

double A_(void){ 

returnC ua * (11(ua)*K0(ub) + 10(ub)*K1(ua)) );} 

double B_(void)<; 

return ( ub'^ub^C 10(ub)*K0(ua)-10(ua)*K0(ub)) / ( 2*k*b_width*depth*m()));} 
double C_(void)<; 

return ( 2*k*a_width*depth*m()*ub*(11(ub)*K1(ua) - 11(ua)*K1(ub))/(ua) );} 
double D_(void)<; 

returnC ub*( I0(ua)*K1(ub)+I1(ub)"^K0(ua)) );} 

^•k ★**★★★***★*★**★**•**★**■★★★★★★★★★★★★***★**★***★*★*★★★★★★*★★★★★★*★★★*★★★★★★*** 


OPTIMIZE Analysis Program 
Longitudinal Fin of Rectangular Profile 


kkkkkkkkkkkkirk'kkk'kirifir'k'k'kkkit'kifiliriririfkkkkititifkkifitifk'kkkitkkkif'kkififir'kirkkkkkitifk'itit'kkkk ★ j 


#include 

<stdio.h> 


#include 

<stdlib.h> 


#include 

<math.h> 


#include 

<conio.h> 


#include 

"units.h" 


#define 

Yes 

1 

#define 

No 

0 

#define 

N 

100 

#define 

s i gma 

5.66961e-8 

double 

length 

, width, depth 


100 




float Tbase, heat, PN, eta, V, K2, epsilon, k, density, mass, Ap, error; 

int T, L, M, S, H, d, K, i, X, 2, Y; 

void inputO, computeO, resultsC); 

★★*★★★***★★★★***★★*★★***★★★★****★**★★**★★**★★★*★★*★★★★*★★*★★★★★****★*★★ 


MAIN PROGRAM 


int main(){ 

FILE *data_in; 

data_in = fopenC'units.ext","r"); 

fscanf(data_in,"%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,&d,&K,&X,&2,&Y,&error); 

fclose(data_in); 

error = error / 100; 

input(); 

heat = heat / ( 1 - V); 

Ap = (1.6178 * pow(heat/depth , 3)) / ( k * pow(epsi Ion,2) pow(Tbase,9) * pow(sigma,2)) 
width = (1.8648* pow(heat/depth , 2)) / ( k * sigma*epsilon*pow(Tbase,5)); 

computeO; 
results(); 
ex i t ( 0 ) ; 
return 0;> 

j-k ■k-k-k-kkit-k-k'kit-kk-k-kit-k-k-k'kk-k-k-k'kit-k-k-k-kk-k'k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k'k'kitk-k-k-k-kitk-k-k'kit-k-kickick'k'kifk 

COMPUTE ROUTINE 

■k-k-kk-k-k-k-kk-kk-k-k-k-kk-k-ttk-k-k-kk-kk-itk-ttkk-kkk'kk-kk-tfk-kkk-k-k-ifk-tfk-k-k-k-k-tfk-k-kk-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k -k j 


void compute(void)<; 

double A_(), B_(), C_(), D_(), m(), hr(), Yb(); 

double inc, overall_diff = 100, Tmax = 0, Tbv^min = 100000, data[N+1] [6]; 

double A1, B1, Cl, D1, Ttip, Tave, temp, qa, segmented!ff; 
int overflow, possible_overflow = No, number_of_tries=0; 

length = Ap / width; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while (fabs(overall_diff) > ( Tbase * error ) J| i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * width * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == NoX 

segment_diff = 10; 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No)( 

A = A_(); B = B_(); C = C^(); D = 0_(); 
temp =A*Ta+B* qa; 
qb = C * Ta + 0 * qa; 
segment_diff = ( temp - Tb) / 2; 
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Tb = temp + segmented!ff; 

if (Tb > Tbase + 50 11 Tb < OX 

overflow = Yes; possible^overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Trip; 
inc = (Tbv_min - Tmax) / 2 ; 

Ttip = fabs(Tbv_min • inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No){ 

data[i][0] = qb; data[i][1] = Tb; data[i][2] = A; 
data[i][3] = B; data[i][4] = C; data[i][5] = 0; 

Ta = Tb; qa = qb; ++i;}} 

if (number_of_tries > 29){ 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(0);> 

if (overflow == No)C 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb * Tbase) / 2; 

if ( possible_overflow == Yes)<; 
inc = inc / 2; 

if (overall_diff >0) Ttip = Ttip- inc; 
else” Ttip = Ttip + inc;) 

else 

Ttip = Ttip - overall_diff; 
number_of_tries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,TTip),w_T(T,Tb)) 
if (Ttip < Tave) Ttip = Tave;}} 

for (i = N - 1; i >= 1; i--){ 

A1 = A * dataCi] [2] + B * datali] [4]; B1 = A * dataCi] [3] + B * dataCi] [5] ; 

Cl = C * datali] [2] + D datali] [4]; 01 = C * datali] 13] + D * datali] [5]; 

A = A1; B = B1; C=C1; 0=01;} 

Ta = Ttip; 

eta = qb / (2 * sigma * epsilon * length * depth * pow(Tbase,4)); 

PN = length * sqrt(sigma * epsilon * pow(Tbase,3) / (k * width)); 
mass = width * length * depth * density * 1000; 
return;} 

j-k ★★★★★★★★★★★*★★★★★★★★★★★★★★★★★★*★*★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

OUTPUT PROCEOURE 

★ A************************************************************************ -k! 


void results(void)<; 


clrscr(); 

cprintf("\a OPT 1 MUM PARAMETERS"); 

cprintf("\r\n THE LONGITUOINAL FIN OF RECTANGULAR PROFILE"); 

cprintf("\r\n ================================s===:==============================\r\n"); 

cprintf(" Oensity = %4.2f kg / dm^3\r\n",density); 

cprintf(" Fin Oepth = "); p_L(L, depth); 

cprintfC Total Mass = "); p_M(M, mass); 

cprintfC Fin Height = "); p_L(L, length); 

cprintf(" Fin Width = "); p_L(d, width); 

cprintfC Emissivity = %4.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X, Ap); 

cprintfC Profile Number = %4.3f\r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4.3f\r\n",eta); 

cprintfC Base Temperature = "); p_T(T,Tb); 
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cprintfC Edge Temperature = *'); p_T(T,Ta); 

cprintfC" Thermal Conductivity = *'); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC*' ====rr======rr============================================r=======\r\n"); 

cprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = %4.4f B = ",A);pZ(2,B); 

cprintfC C = “); p_Y(Y,C); 

cprintfC D = %4.4f", D); 

getchO; 

return;> 


J-k ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★I'r***************************************** 

INPUT PROCEDURE 


*★**★*★★★★*★★★★***★*★*★★**★*★★*★*★**★★*★★**★*******★*★****★★*★*★★★★★★★★★★★ ★j 


void input(void)( 

clrscrC); _setcursortype(_NORMALCURSOR); 


cprintfC\r\n THE LONGITUDINAL FIN OF RECTANGULAR PROFILE\r\n"); 

cpr intf C\r\n ===================rr=rrrrr==r==r==rrr=r=rr====r==r=rr=r==========\r\n") 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Heat Dissipation "); heat = i_H(H); 

cprintfC Emissivity = "); scanfC%f", &epsilon); 

cprintfC Density ( kg / dm"3 ) = "); scanfC%f", ^density); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanfC%f", &V); 

cprintfC\r\n rrr=====r ================================ ====r========r===r===r ===»«); 

if ( epsilon > 1.0 j| epsilon <= 0)( 


cprintfC\r\n\r\nFATAL ERROR; Emissivity is not between 0 and 1"); 
getchO; exit(O);} 

if ( V < 0 1 I V > 0.8){ 

cprintfC\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 
getchO; exit(0);> 

_setcursortype(_N0CURS0R); 
return;> 

j-k kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

GLOBAL FUNCTIONS 

kkkkkkkk-kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk k j 


double hr(void){ 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 

3) - (K2 / 

C2 * Tav)));} 

double m(void){ 

returnC sqrt( ( 2.0 * hr()) / ( k 

* width) ) 

);> 

double Yb(void)( 

returnC sqrt( 2.0 * hr() * k * width) * depth 

);> 

double A_(void)( 

returnC coshC mC) * length / N ) 

);> 


double B_Cvoid){ 

return C sinhC mC) * length / N ) 

/ Ybo );> 


double C_Cvoid)( 

return C sinhC mC) * length / N ) 

* YbO );> 


double D_Cvoid){ 
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returnC cosh( m() * length / N ) 


);> 
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jtt *********i(*i(*************tc-ktfk-ktttc**it**-k-k*iiiit^ifk*ii 

OPTIMIZE Analysis Program 
Longitudinal Fin of Trapezoidal Profile 

*******************»*******************)k*******)k:m 


****************** 


***★★★**★***★**★*★<» 


#include 

<stdio.h> 


#include 

<stdlib.h> 


#include 

<math.h> 


#include 

<conio.h> 


#include 

"units.h" 


#include 

"bessel.h" 


#define 

Yes 

1 

#define 

No 

0 

#define 

N 

100 

#define 

s i gma 

5.66961e-8 

char 

file name[20]; 

double 

length 

, edge_width 

double 

A, B, 

C, D, Ta, a. 

float 

K2, taper_ratio, V 

int 

T, L, 

M, S, H, d, 1 

void 

inputO, computeO 

FILE 

*data_ 

in, *out; 


*********4 


fc******************************!*:*************** 


MAIN PROGRAM 


★ ★★★★★^^★★★★★★★★★★★★★♦★★★★★★★★★★★★★★★ik^t***************************** 


int mainCX 

data_in = fopen("uni ts.ext*', "r"); 

fscanf (data_in/'%d %d 7od 7pd %d %d %d %d %d %d %f",&T ,&L,&M,&S,&H,&d,&<,&X,&Z,&Y,&error); 
fclose(data_in); 
error = error / 100; 

inputC); 

heat = heat / (1 - V); 

Constant_3 = 0.3881*powCtaper_ratio,2) - 1.0324*taper_ratio + 2.7522; 

Ap = ((0.5147 * taper_ratio + 1,1382) * pow(heat/depth , 3)) / ( k * pow(epsilon,2) * 

pow(Tbase,9) * pow(sigma,2)); 

base_width = (Constant_3 * pow(heat/depth , 2)) / ( k * sigma*epsilon*pow(Tbase,5)); 
edge_width = taper_ratio * base_width; 

compute(); 
results(); 
exit(O); 
return 0;> 

y* **********★**★★*********★****★***********★★*******★*★********★★★★**★**★ 

COMPUTE ROUTINE 

************************************************************************** *y 

void compute(void){ 

double A_(), B_(), C_(), D_(), m(), hr(), u(double x); 

double length_inc, bT, db; 

double inc, Tave, overall^diff = 100, Tmax = 0, Tbv_min = 100000, data [N+1] [6]; 

double Al, B1, Cl, D1, Ttip, temp, qa, segment_diff; 

int overflow, possible_overflow = No, number_of_tries=0; 
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length = 2 * Ap / ( base_width + edge_width ); 
b1 = length / ( 1 - taper^ratio ); 
fin_taper = atan2( base_width, 2*b1 ); 
db = length / N; 

K2 = 2 * epsilon * sigma * V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=1; 

while (fabs(overall_diff) > ( Tbase * error) 1] i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = b1 - length; 

while (i <= N && overflow == No)<; 
segment_diff = 10; 

a_width = edge_width * a / ( b1 - length); 
b_width = edge_width * ( a + db ) / ( b1 - length); 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No)<; 

ua = u(a) ; ub = u(a + db); 

A = A_(); B = B_(); C = C_(); D = D_(); 
temp =A*Ta+B*qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase + 50 || Tb < 0)C 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No)<; 

data[i][0] = qb; data[i3[1] = Tb; data[i][2] = A; 
data[i][3] = B; data[i][4] = C; data[i][5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}} 

if (number_of_tries > 29){ 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(O);} 

if (overflow == No)<; 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes)<; 

inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overall_diff; 
number_of_t ries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)) 

if (Ttip < Tave) Ttip = Tave;}} 
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for (i = N - 1; i >= 1; i-‘)( 

AT = A * data[i][2] + B * data[i][A]; B1 = A * dataCi]C3] + B * data[i][5]; 

Cl = C * dataCi] [2] D * dataCi] CA]; Dl = C * data[i] C3] + D datafi] [5]; 

A = A1; B = Bl; C=C1; D=D1;) 

if (S == DC 

out = fopen(file_name, "w"); 
for (i = 1; i <= N; i++) 

fprintfCout, "%e, %e, %e\n", w_L(d, i * length / N), w_T(T,dataCi][1]), w_H(H, 

datali] [0] ) ); 

fclose(out);> 

Ta = Ttip; 

Ap = .5 * length * ( edge_width + base_width ); 
mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow(length, 2) * pow(Tbase,3) / (k* base_width); 
eta = qb / ( 2 sigma * epsilon * length * depth * pow(Tbase,A)); 
return;) 

OUTPUT PROCEDURE 


void results(void){ 


cIrscr(); 

cprintf("\a\r\n OPTIMUM PARAMETERS"); 

cprintf("\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE"); 

cprintf("\r\n ==================================================================\r\n"); 

cprintf(" Density = %A.2f kg / dm*3\r\n",density); 

cprintfC" Fin Depth = "); p_L(L, depth); 

cprintfC Total Mass = "); p_M(M, mass); 

cprintfC Fin Length = "); p_L(L, length); 

cprintfC Emissivity = %A.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X,Ap); 

cprintfC Profile Number = %A.3f\r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %A.3f\r\n",eta); 

cprintfC Fin Edge Width = "); p_L(d, edge_width); 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC Base Temperature = "); p_T(T,Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %A.2f\r\n", V); 

cprintfC = = = = = = = = r = r = = = s: = = rr = rrr = = r = r = = = = = = = = = = = r = = = = = = = = = = = = = = = = = = = = = = = = = = \ r\n" ) ; 

CprintfC Transmission Parameter Matrix \r\n\r\n"); 

cprintfC A = U.Af B = ", A); p Z(Z,B); 

cprintfC C = "); p_Y(Y,C); 

cprintfC D = %A.Af", D); 

getchO; 

return;} 


j* *-k •*(-k *-k-k-k-k -k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k it-k **-k-k-k-k-k •*(-k it ick-ick-ick-k-kickick-k-k-k-k-k •*(■*(•»(•»( it-kick 


INPUT PROCEDURE 

*★★★★**★★★**********★*****★*******★*★*******************★★★★*★★★★★*★**★*■★★ k j 


void input(void)( 

clrscr(); _setcursortype(_NORMALCURSOR); 


cprintfC\r\n THE LONGITUDINAL FIN OF TRAPEZOIDAL PROFILE\r\n"); 

cprintfC\r\n ========r=========================================================\r\n"); 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Heat Dissipation "); heat = i_H(H); 

cprintfC Taper Ratio = "); scanfC%f",&taper_ratio); 

cprintfC Emissivity = "); scanfC%f", Sepsilon); 
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cprintfC" Density ( kg / dm"3 ) = "); scanf("%f", &density); 

cprintfC" Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanfC%f", &V); 

cprintfC\n==================================================================\r\n"); 


if ( epsilon > 1.0 ]] epsilon <= OX 

cprintf(“XrXnXrXnFATAL ERROR: Emissivity is not between 0 and 1"); 
getchO; exit(0);> 

if ( V < 0 1j V > 0.8){ 

cprintfCXrXnXrXnFATAL ERROR: Environmental Parameter is not between 0 and 0.8"); 
getchO; exit(0);> 

if ( taper_ratio > .99 ]] taper_ratio <= 0)<; 

cprintfCXrXnXrXnFATAL ERROR: Taper Ratio is not between 0 and 1.0"); 
getchO; exit(0);> 

_setcursortype(_NOCURSOR); 

return;} 

***********★★*****★★★**★*★******★**★★***★★★*★****★**★★★★*★**★★★*★★*★★*★ 

GLOBAL FUNCTIONS 

************************************************************************** *j 


double hr(void){ 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));> 
double m(void){ 

returnC hr() / (k * sin(fin_taper)));> 

double uCdouble x){ 

returnC 2 * sqrtCx * m() ) );> 

double A_(void)( 

returnC ua * C11Cua)*K0Cub) + 10Cub)*KlCua)) );> 

double B_Cvoid){ 

return C ub*ub*CI0Cub)*K0Cua)-I0Cua)*K0Cub)) / C 2*k*b_width*depth*mC)));> 
double C_Cvoid)<; 

return C 2*k*a_width*depth*mC)*ub*C11Cub)*KlCua) - 11Cua)*K1Cub))/Cua) );> 
double D_Cvoid){ 

returnC ub*C I0Cua)*KlCub)+I1Cub)*K0Cua)) );> 
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j-k ** 


************************** 


************************* 


********* 


OPTIMIZE Analysis Program 
Longitudinal Fin of Triangular Profile 

**************************************************************************** kj 


#inc1ude 

<stdio.h> 


#include 

<stdlib.h> 


#include 

<math.h> 


#include 

<conio.h> 


#include 

"units.h" 


#include 

"bessel.h" 


#define 

Yes 1 


#define 

No 0 


#define 

N 99 


#define 

sigma 5.66961e-8 


char 

file_name [20]; 


double 

length, edge width. 

base_width 

double 

A, B, C, D, Ta, Tb, 

qb=0, ua, i 

float 

V, K2, epsilon, k. 

Tbase, dens 

int 

T, L, M, S, H, d, K, i, X, 2, 

void 

inputO, results(), 

compute(); 

FILE 

*data_in, *out; 



PN, heat; 


/* *************************** 


**************************************** 


MAIN PROGRAM 


************************ 


***************************************** */ 


int main(){ 

data_in = fopenC'units.ext","r"); 

fscanf(data_in/'%d %d %d %d %d %d %d %d %d %d %f",&T,&L,&M,&S,&H,terror); 

fclose(data_in); 

error = error / 100; 

inputC); 

heat = heat / (1 - V); 

Ap = (1.1234 * pow(heat/depth , 3)) / ( k * pow(sigma*epsiIon,2) * pow(Tbase,9)); 

base_width = (2.4852 * pow(heat/depth , 2)) / ( k * sigma*epsilon*pow(Tbase,5)); 

compute(); 
results(); 
ex i t ( 0 ) ; 
return 0;) 

j* *********************************************************************** 


COMPUTE ROUTINE 

************************************************************************** *y 


void computeOC 


double 

double 

double 

double 

int 


A_(), B_(), C_(), D_(), m(), hr(), u(doubte x); 
length_inc, bl, db; 

inc, Tave, overall^diff = 100, Tmax=0, Tbv^min = 100000, data[N+1][6]; 
Al, Bl, Cl, Dl, Ttip, temp, qa, segment_diff; 
overflow, possible_overflow = No, number_of_tries = 0; 


length = 2 * Ap / base_width; 
edge_width = base_width / 100; 
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db = length / ( N 1 ); 

bl = length; 

length = length - db; 

fin_taper = atan2( base_width, 2^b1 ); 

K2 = 2 * epsilon * sigma ^ V * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 ^ Tbase * sqrt( 2 ^ V ) + ( ( Tbase - (Tbase sqrt( 2 * V )))/2) 

else 

Ttip = Tave + 0.0001; 
inc = .5 * Ttip; i=l; 

while (fabsCoverall_diff) > ( Tbase * error) }j i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segment_diff = 10; 

while (fabs(segment_di f f) > ( Tb ^ error) / N && overflow == No)<: 

ub = u(db); 

temp = Ta / lO(ub); 

C = 2 * k * edge_width * depth * pow(m(),2) ^ Il(ub) / (ub * lO(ub)) 
qb = C ^ temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp **■ segment_di f f ; 

if (Tb > Tbase 50 | 1 Tb < 0)<; 

overflow = Yes; possible_overfIow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No)<; 

data[i] [0] = qb; data[i][l] = Tb; data[i][2] = 1; 
data[i][3] = 0; data[i][4] = C; dataCi] [5] = 1; 

Ta = Tb; qa = qb; ++i;} 


a = bl - length; 

while (i <= N && overflow == No)( 
segment_diff = 10; 

a_width = edge_width ^ a / ( bl - length); 
b_width = edge_width ’^ ( a + db ) / ( bl - length); 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No){ 

ua = u(a) ; ub = u(a + db); 

A = A_(); B = B^O; C = C_(); 0 = D_(); 
temp = A^Ta-*-B^ qa; 
qb =C*Ta + D* qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase 50 j 1 Tb < 0){ 

overflow = Yes; possible_overfIow = Yes; 
if (Ttip < Tbv__min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < TaveT Ttip = Tave;}} 

if (overflow == No){ 

data[i][03 = qb; data[i][l] = Tb; dataCi] [2] = A; 
data[i][3] = B; dataCi] [4] = C; dataCi] [5] = D; 

Ta = Tb; qa = qb; a = a + db; ++i;}} 
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if (number_of_tries > 29)( 

clrscrO; cprintf("Fatal Error: Unable to Compute"); 
getchO; exit(O);) 

if (overflow == NoX 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == YesTC 
inc = inc / 2; 

if (overal l__di f f >0) Ttip = Ttip* inc; 
else Ttip = Ttip + inc;) 

else 

Ttip = Ttip * overall_diff; 
number_of_t ries++; 
gotoxy(27,24); 

cprintfC'N = %d Ta = %4.1f Tb = %4.1f ",number_of_tries,w_T(T,Ttip),w_T(T,Tb)); 
if (Ttip < Tave) Ttip = Tave;)) 

for (i = N - 1; i >= 1; i‘-)C 

A1 = A * data[i][2] + B * data[i][4]; B1 = A * data[i][3] + B * data[i][53; 

Cl = C * data[i][2] + D * data[i][4]; D1 = C * data[i][3] + D * data[i][5]; 

A = AT; B = B1; C=C1; D=D1;} 

length = length + db; 

if (S == 1)( 

out = fopen(file_name, "w"); 
for (i = 1; i <= N; i++) 

fprintf(out, "%e, %e, %e\n", w_L(d, i * length / N), w_T(T,data [i3 [13), w_H(H, 

data[i3 [03)); 

fclose(out);} 


Ta = Ttip; 

mass = Ap * depth * density * 1000; 

PN = 2 * sigma * epsilon * pow(length,2) * pow(Tbase,3) / (k*base_width); 
eta = qb / ( 2 * sigma * epsilon * length * depth * pow(Tbase,4)); 
return;) 

/* *****************************************^**************************** 


OUTPUT PROCEDURE 


void results(void)(; 


c I rscr( ); 

cprintf("\a OPTIMUM PARAMETERS"); 

cprintf("\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROFILE"); 

cprintf("\r\n ==================================r===============================\r\n"); 

cprintfC Density = %4.2f kg / dm^3\r\n",density); 

cprintfC Fin Depth = "); p_L(L, depth); 

cprintfC Total Mass = "); p_M(M, mass); 

cprintfC Fin Length = "); p_L(L, length); 

cprintfC Emissivity = %4.2f\r\n", epsilon); 

cprintfC Profile Area = "); p_A(X,Ap); 

cprintfC Profile Number = %4.3f\r\n",PN); 

cprintfC Base Heat Flow = "); p_H(H, qb); 

cprintfC Fin Efficiency = %4.3f\r\n",eta); 

cprintfC Fin Base Width = "); p_L(d, base_width); 

cprintfC Base Temperature = "); p_T(T,Tb); 

cprintfC Edge Temperature = "); p_T(T,Ta); 

cprintfC Thermal Conductivity = "); p_K(K, k); 

cprintfC Environmental Parameter = %4.2f\r\n", V); 

cprintfC ==================================================================\r\n"); 
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cprintf(" 
cprintf(" 
cprintf(" 
cprintf(" 
getchC); 
return;} 


Transmission Parameter Matrix \r\n\r\n“); 
A = %4.4f B = ", A); p_2(2,B); 

C = "); p_Y(Y,C); 

D = %A.4f", D); 


j-k ■k-k-k-k-k-k-kkk'k-k-kk-k-k-k-k-k-k'k-k-k-k-kk-k-k-k-k'k'k'k'k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k-k-k-k-k-kk-k-k-k'k'k'k-k-k-k-k-k-kk-k-kk 


INPUT PROCEDURE 


■k-kkkk-k-kkkkkkk-k-k-k-k-k-kkkkk-k-kk-k-k-k-k-k-k-k-k-k-k-k-k-k-k-kk-kkk-k-k-k-kk-k-k-k-k-k-k-k-k-kkkk-kkk-k-k-tck-k-kkk-k k j 


void input(void){ 

clrscr();_setcursortype(_NORMALCURSOR); 


cprintf("\r\n THE LONGITUDINAL FIN OF TRIANGULAR PROFILE"); 

cprintf("\r\n ===========================================================:======r\r\n") 

cprintfC Fin Depth "); depth = i_L(L); 

cprintfC Heat Dissipation "); heat = i_H(H); 

cprintfC Emissivity = "); scanfC%f", &epsilon); 

cprintfC Density ( kg / dm^3 ) = "); scanf("%f", Sdensity); 

cprintfC Base Temperature "); Tbase = i_T(T); 

cprintfC Thermal Conductivity "); k = i_K(K); 

cprintfC Environmental Parameter (0 to .8) = "); scanfC%f", SV); 

cprintfC\n==================================================================\r\n"); 

if ( epsilon > 1.0 || epsilon <= 0){ 


cprintfC\r\n\r\nFATAL ERROR: Emissivity is not between 0 and 1"); 
getchO; exit(0);> 

if ( V < 0 11 V > 0.8){ 

cprintf("\r\n\r\nFATAL ERROR: Environmental Parameter is not between 0 and 0,8"); 
getchO; exit(O);} 

_setcursortype(_NOCURSOR); 
return;> 

j k kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk-kkkkkkk-kk-kkkkkkkkkkkkkkkkkkkkkkkkkkkk 

GLOBAL FUNCTIONS 


*****★**★******★********************★***★★★*****★★■***★★*★*★***★***★***★**★ k j 


double hr(void)( 

float Tav = (Tb + Ta) / 2; 

returnCsigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)));} 
double m(void)( 

returnC hr() / (k * sin(fin_taper)));} 

double uCdouble x){ 

returnC 2 * sqrtCx * mC) ) );> 

double A_Cvoid){ 

returnC ua * C11(ua)*K0Cub) + I0Cub)*KlCua)) );> 

double B_Cvoid){ 

return C ub*ub*CI0Cub)*K0Cua)-I0Cua)*K0Cub)) / C 2*k*b_width*depth*mC)));> 
double C_Cvoid){ 

return C 2*k*a__width*depth*mC)*ub*C11Cub)*KlCua) - 11Cua)*KlCub))/Cua) );} 
double D_Cvoid)<; 

returnC ub*C I0Cua)*K1Cub)+IlCub)*K0Cua)) );} 
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j-k *************************************************************************** 


UNIT CONVERSION Program Header File 


*****************************************«*******««*******4r*********ir******« k j 

^include <stdio.h> 

JS^include <conio.h> 

^include <string.h> 

^include <stdlib.h> 

^include "units.h" 

float i_T(int T 
float V; 
switch (T)( 

case 0: 
case 1: 
case 2; 
case 3: 
returnCV);} 

void p_T(int T, 
switch (T){ 

case 0: 
case 1: 
case 2: 
case 3: 
return;> 

float w_T(int T, 
switch Tt)( 

case 0: 
case 1: 
case 2: 
case 3: 
return(V);} 

void p_Z(int 1, float V)( 
switch (Z)( 

case 0: cprintf("%4.4f deg C / watt\r\n", V); break; 

case 1: V = 0.1760831 * V ; cprfntf("%4.34 deg F hr / BTU\r\n", V);} 

return;} 

float w_Z(int Z, float V){ 
switch (Z)( 

case 0: break; 

case 1: V = 0.1760831 * V;} 

return(V);} 


float V)( 
break; 

V = V - 273.15; break; 

V = 9 * ( V - 255.222222) / 5; break; 

V = 9 * V / 5;} 


)( 


cprintf("(deg K) = "); scanf("%f", &V); break; 
cprintf("(deg C) = "); scanf("%f", &V); V = V + 273.15; break; 
cprintf("(deg F) = "); scanf("%f", &V); V = 5 ^ ( V + 459.4 ) / 9; 
cprintf ("(deg R) = "); scanf("7of", &V); V = 5 * V / 9;} 


float V)( 


cprintf("%4.2f deg K\r\n‘', V); break; 

V = V ' 273.15; cprintf("%4.2f deg C\r\n", V); break; 

V = 9 * ( V - 255.222222) / 5; cprintf("%4.2f deg F\r\n", V); 

V = 9 * V / 5; cprintf("%4.2f deg R\r\n", V);} 


break; 


void p_Y(int Y, float V)<; 
switch (Y)< 

case 0: cprintf("%4.4f watt / deg C", V); break; 

case 1: V = V / 0.1760831 ; cprintf("7o4.4f BTU / hr deg F", V);} 

return;} 

float w_Y(int Y, float V)( 
switch (Y){ 

case 0: break; 

case 1: V = V / 0.1760831;) 

return(V);} 

float i_L(int L){ 
float V; 
switch U){ 


case 

0: 

cprintf(" 

(m) 

= "); 

scanf("7of", 

&V); 

break; 

case 

1: 

cprintf(" 

(cm) 

= "); 

scanf("%f", 

&V); 

V = 

V / 100; break; 

case 

2: 

cprintf(" 

(ft) 

= "); 

scanf("%f", 

&V); 

V = 

0.3048 * V; break; 

case 

3: 

cprintf(" 

(in) 

= "); 

scanf("%f", 

&V); 

V = 

0.0254 * V; break; 

case 

4: 

cprintf(" 

(mm) 

= "); 

scanf("%f", 

&V); 

V = 

V / 1000;} 


break 
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return(V);} 


void p_L(int L, float V)<; 
switch (L)< 

case 0: cprMitf("%4,2f m\r\n", V); break; 

case 1: V = 100 * V; cprintf("%4.2f cm\r\n", V); break; 

case 2: V = V / 0.3048; cprintf("%4.2f ft\r\n", V); break; 

case 3: V = V / 0.0254; cprintf("%4.2f in\r\n", V); break; 

case 4: V = 1000 * V; cprintf("%4.2f mm\r\n", V);} 

return;} 


float w_L(int L, 
switch (L){ 

case 0: 
case 1: 
case 2: 
case 3; 
case 4: 
return(V);} 


float V)<; 
break; 

V = 100 * V; break; 

V = V / 0.3048; break; 

V = V / 0.0254; break; 

V = 1000 * V;} 


float i_A(int A){ 
float V; 
switch (A)<; 


case 

0: 

cprintf( 

*' (m"2) 

= "); 

scanf("%f ”, 

&V); break; 



case 

1: 

cprintf( 

"(cm"2) 

= ”); 

scanf C%f*', 

&V); V = V 

/ 

10000; break; 

case 

2: 

cprintf( 

"(ft^2) 

= "); 

scanf(“%f". 

&V); V = V 

/ 

10.76; break; 

case 

3: 

cprintf( 

"(in"2) 


scanf("%f", 

&V); V = V 

/ 

1550; break; 

case 

4: 

cprintf( 

"(rT¥n^2) 


scanf("%f", 

&V); V = V 

/ 

1000000;} 


returnCV);) 


void p_A(fnt A, float V)< 
switch (A)<; 

case 0: cprintf("%4.2 m " 2\r\n", V); break; 

case 1: V = 10000 * V; cprintf(*'%4.2f cm " 2\r\n", V); break; 

case 2: V = V * 10.76; cprintf("%4.2f ft ^ 2\r\n'‘, V); break; 

case 3: V = V * 1550; cprintf("%4.2f in ^ 2\r\n", V); break; 

case 4: V = 1000000 * V; cprintf("%4.2f mm ^ 2\r\n“, V);} 

return;} 


float w_A(int A, 
switch (A)<; 

case 0: 
case 1: 
case 2: 
case 3; 
case 4: 
return(V);} 


float V){ 
break; 

V = 10000 * V; break; 

V = V * 10.76; break; 

V = V * 1550; break; 

V = 1000000 * V;} 


float i_^M(int M){ 
float V; 
switch (M){ 

case 0; 
case 1: 
case 2: 
case 3: 
return(V);} 


cprintfC (kg) 
cprintfC (gm) 
cprintf("(Ibm) 
cprintfC (oz) 


"); scanfC%f", 
"); scanfC%f", 
"); scanf("%f", 
"); scanf("%f'‘, 


&V); break; 

&V); V = V / 1000; break; 

&V); V = V / 2.204622; break; 
&V); V = V / 35.27396; break;} 


VOid p_M (int M, 

switch (M)<; 

case 0: 
case 1: 
case 2: 
case 3: 
return;} 


float V){ 

cprintf("%4.2f kg\r\n", V); break; 

V = 1000 * V; cprintf("%4.2f grams\r\n", V); break; 

V = 2.204622 * V; cprintfC%4.2f lbm\r\n", V); break; 

V = 35.27396 * V; cprintfC7o4.2f 02 \r\n", V);} 


float w_M(int M, 
switch Tm){ 

case 0: 
case 1: 


float V)<; 


break; 

V = 1000 * V; break; 
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case 2: V = 2.204622 * V; break; 

case 3: V = 35.27396 * V;} 

return(V);} 


float i_H(int H){ 
float V; 
switch (H)( 

case 0: 
case 1: 
case 2: 
case 3: 
returnCV);} 


cprintf(“(watts) = 
cprintfC (kw) 
cprintf(“(BTU/sec) 
cprintfC (BTU/hr) 


“); scanfC%f“, &\/); break; 

= “); scanfC%f“, &V); V = 1000 * V; break; 

= “); scanfC%f“, &V); V = V * 1054.3502; break; 

= “); scanf(“%f“, &V); V = V / 3.414425;} 


void p_H(int H, 
switch (H)( 

case 0: 
case 1; 
case 2: 
case 3: 
return;} 


float V){ 

cprintf(“%4.2f watts\r\n“, V); break; 

V = V / 1000; cprintf(“%4.2f kw\r\n'', V); break; 

V = V / 1054.3502; cprintf(“%4.2f BTU / sec\r\n“, \/);break; 

V = 3.414425 * V; cprintf(“%4.2f BTU / hr\r\n“. V);} 


float w_H(int H, 
switch (H){ 

case 0: 
case 1: 
case 2: 
case 3: 
return(\/);} 


float VX 
break; 

V = V / 1000; break; 

V = V / 1054.3502; break; 

V = 3.414425 * V;} 


float i_K(int K)( 
float V; 
switch (K)( 

case 0: cprintf (“(watts / m deg K) = “); scanf(“%f“, &\/); break; 

case 1: cprintf (“(BTU / ft hr deg F) = “); scanf(“%f“, &\/); V = V / 1.731;} 

return(\/);} 


void p_K(int K, float V){ 
switch (KX 

case 0; cprintf(“%4.2f watts / m deg K\r\n“, V); break; 

case 1: V = 1.731 * V; cprintf(“%4.2f BTU / ft hr deg F\r\n“, V);} 

return;} 
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j* ************************************************************************** 


BESSEL FUNCTION Program Header File’^ 


*************************************************************************** 

^include <stdio.h> 

#include <math.h> 

#include "bessel.h" 

double K0(float x)( 

double y, z; 

double Pl=-0.57721566, P2=0.4227842, P3=0.23069756, P4=0.0348859, P5=0.00262698, 

P6=0.0001075, P7=0.0000074; 

double 01=1.25331414, Q2=-0.07832358, 03=0.02189568, Q4=-0.01062446, 05=0.00587872, 

06=-0.00251546; 
double 07=0.00053208; 

if ( X < 2 ){ 

y = ( X * X ) / 4; 

z = (-log(x/2)*I0(x))+(PI+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7))))));) 

elseC 

y = 2 / x; 

z = (exp(-x)/sqrt(x))*(0l+y*(02+y*(03+y*(04+y*(05+y*(06+y*07))))));} 
returnC 2 );) 

double KKfloat x){ 

double y, z; 

double P1=1.0, P2=0.15443144, P3=-.67278579, P4=-0.18156897, P5=-0.01919402, P6=-0.00110404, 

P7=-0.00004686; 

double 01=1.25331414, 02=0.23498619, O3=-0.03655620, 04=0.01504268, 05=-0.00780353, 

06=0.00325614; 

double Q7=-0.00068245; 

if ( X < 2 ){ 

y = ( X * X ) / 4; 

z = (log(x/2)*11(x))+(1/x)*(Pl+y*(P2+y*(P3+y*(P4+y*(P5+y*(P6+y*P7))))));} 

else( 

y = 2 / x; 

z = (exp(-x)/sqrt(x))*(Ql+y*(Q2+y*(Q3+y*(Q4+y*(Q5+y*CQ6+y*Q7))))));} 
return(z);} 


double I0(float x)<. 
double y, z, ax; 

double P1=1, P2=3.5156229, P3=3.0899424, P4=1.2067492, P5=0.2659732, P6=0.0360768, 

P7=0.0045813; 

double 01=0.39894228, 02=0.01328592, 03=0.00225319, 04=-0.00157565, 05=0.00916281, 

Q6=-0.02057706; 

double 07=0.02635537, O8=-0.01647633, 09=0.00392377; 

if (fabs(x) < 3.75 )<; 

y = ( X / 3.75 ) * ( X / 3.75 ); 
z = Pl+y*(P2'*-y*(P3+y*(P4'^y*(P5 + y*(P6+y*P7)))));) 

elseL 

ax = fabs(x); 
y = 3.75 / ax; 

z = (exp(ax)/sqrt(ax))*(0l+y*C02+y*(03+y*(04+y*(05+y*C06+y*(07+y*(08+y*09))))))));> 
return(z);} 

double IKfloat x){ 

double y, z, ax; 

double P1=0.5, P2=0.87890594, P3=0.51498869, P4=0.15084934, P5=0.02658733, P6=0.00301532, 

P7=0.00032411; 

double 01=0.39894228, 02=-0.03988024, o3=-0.00362018, 04=0.00163801, Q5=-0.01031555, 
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06=0.02282967; 

double Q7=-0.02895312, 08=0.01787654, o9=-0.00420059; 

if (fabs(x) < 3.75 )( 

y= ( x/ 3.75 ) * ( x/ 3.75 ); 

2 = x*(Pl + y*(P2+y*(P3 + y*(P4-*-y*(P5+y*(P6+y*P7))))));) 

else( 

ax = fabs(x); 
y = 3.75 / ax; 

2 = (exp(ax)/sqrt(ax))*(0l+y*(02+y*(03+y*(04+y*(05+y*(06+y’^(Q7+y*(Q8+y*Q9))))))));> 
return( 2 );} 
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APPENDIX B COMPUTER CODE TO DERIVE OPTIMUM GEOMETRY 


I* ****★★★★★★★★**★**★★**********★*★**★★***★★★★★★********★***★★★★★★*★★★*★★★**★★ 

PROFILE AREA Analysis Program To Determine The Optimum Geometry 
Longitudinal Fin of Rectangular Profile 


#include <time.h> 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <conio.h> 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 100 

sigma 5.66961e-8 


double 
double 
float 
i nt 
void 


length, width_base, max_length; 

k, B, C, D, Ta, Tb, qb=0, Heat, qb_max, db, Continue; 

V, K2, epsilon, k, Tbase, Ap, taper_ratio, depth, error=0.01; 
i, j, number_of_passes, failures; 
compute(); 


j-tt ***★*★★★★*******★★**★**★*********************************************** 


MAIN PROGRAM 


★ ★*★**★*★★**★★******★★★**★★******★*******★*★*★****★**★★*★★★★★★**★★★★**★★★★ it j 

int mainOG 
float dq; 

FILE *out; 
randomi2e(); 
error = error / 100; 
cIrscr() ; 

for (j=1; j<102; j++){ 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = .01 + epsilon / 100; 

depth = (float) (20 + random(200)); 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(500); 

Heat = (float) (25 + random(IOOO)); 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 100; 

V = 0; 

taper_ratio = (float) (random(98)); 
taper_ratio = 0.01 + taper_ratio / 100; 
taper_ratio = 1; 


Continue = Yes; 
qb max = * 1; 
qb"= 0; 
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dq = 0; 

Ap = 1.635 * pow(Heat,3) / ( k * pow(sigma,2) * pow(epsiIon,2) * pow(Tbase,9)); 

length = .87 * Heat / ( sigma * epsilon * pow(Tbase,4)); 

db = length / 20016; 

length = length - db; 

width_base = Ap / length; 

if (depth < 10*width_base) Continue = No; 
if (length < 10*width_base) Continue = No; 

while (Continue == Yes){ 

length = length + db; 
width_base = Ap / length; 
number_of_passes+-*-; 
failures = 0; 

compute(); 

while (qb == 0 && Continue == Yes){ 
failures++; 

length = length + db / 2; 
width_base = Ap / length; 
compute(); 
if (failures == 3)C 
Continue = No; 
number_of_passes = 1;}} 


dq = qb - dq; 

cprintf("\r\nwidth_base = %4.3f length = %4.6f heat = %4.6f dq = %4.6f*', 
width_base*1000, length, qb, dq); 
dq = qb; 


if (qb_max < qb ){ 

max_length = length; 
qb_max = qb;} 

else 

Continue = No;} 

if (number_of_passes > 2)C 
length = max_length; 
width_base = Ap / length; 
compute(); 

cprintf ("\r\nWri ting Rec # %d d = %4.3f b = %4.6f qb = %4.3f Tb = %4.1f Ap = %4.1f'‘, 
j, width base, length, qb, Tbase, Ap); 
out = fopen(‘'OPT_RECT.MAT",‘'a‘’); 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e\r\n", j, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper^ratio, width_base, width^base ) 
fclose(out);} 

elseC 

cprintf("\r\nIgnoring This Set Of Parameters"); 

return 0;} 

y* *★★★★★**★★★★★****★*★***★★★★***★*★***★**★*★★★★****★*★****★★★★****★***★★★ 

COMPUTE ROUTINE 


*************************************************************************«: 


void compute(void)<; 

double Tav, inc, overall_diff = 100, Tmax = 0, Tbv_min = 100000; 

double A, B, C, D, m, hr, Yo, Ttip, temp, qa, Tave, segment^diff; 

int overflow, possible_overflow = No, number_of_tries=0; 

width_base = Ap / length; 
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K2 = 2 * epsilon * sigma * \j * pow(Tbase,4); ’ 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Trip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 

else 

Ttip = .0001 + Tave; 
inc = .5 * Ttip; i = 1; 

while (fabs(overall_diff) > ( Tbase * error) I I i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * width__base * epsilon * pow(Ta, 4) * depth; 

while (i <= N && overflow == No)<; 

segment_diff = 10; 

while (fabs(segment_di f f) > ( Tb error) / N && overflow == No)<; 

Tav = (Tb+ Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
m = sqrt( ( 2.0 * hr) / ( k * width_base) ); 

Yo = sqrt( 2.0 * hr k * width_base) * depth; 

A = cosh( m * length / N ); 

B = sinh( m * length / N ) / Yo; 

C = sinh( m * length / N ) * Yo; 

D = A; 

temp =A*Ta+B* qa; 
qb = C * Ta D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase + 50 1 | Tb < 0)<: 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv__min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No){ 

Ta = Tb; qa = qb; ++i;}} 


if (number_of_tries > 29)<; 

Ta=0;Tb=0;qb=0;return;} 

if (overflow == No)( 

if (Tmax < Ttip) 

if (Tb < Tbase ) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible__overf low == Yes)( 
inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip - inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overall_diff; 
number_of_t ries++; 
if (Ttip < Tave) Ttip = Tave;}} 

return;} 
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j-it *A*****i*r************************»**********************************ifc******* 

PROFILE AREA Analysis Program To Determine The Optimum Geometry 
Longitudinal Fin of Trapezoidal Profile 

★ ★★★★★********************************************************************** -kj 


#include <time.h> 
^include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 

#include <conio.h> 
^include "bessel.h" 


#define 
#define 
#define 
#define 


Yes 1 

No 0 

N 100 

sigma 5.66961e-8 


double max_length, length, depth, edge_width, base_width, Ta, Tb, qb=0, inc, Tave; 

float dq, V, qb_max, db, epsilon, taper_ratio, k, error=0.01, Tbase, K2, Heat, Ap; 

int i. Continue, number__of_passes, j, Exit, failures; 

void computeO; 


I* *********************************************************************** 


MAIN PROGRAM 


************************************************************************** * j 


int main(){ 

FILE *out; 
random!zeC); 
error = error / 100; 
clrscrC); 

for (j=1; j < 102; 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = ,01 + epsilon / 100; 

depth = (float) (20 + random(200)); 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(500); 

Heat = (float) (25 + random(1000)); 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 100; 

V = 0; 

taper_ratio = (float) (random(98)); 
taper_ratio = 0.01 + taper_ratio / 100; 

Continue = Yes; 

Exit = No; 
qb_max = -1; 
dq = 0; 
qb = 0; 

Ap = (.63077 * taper^ratio + 1.31723) * pow(Heat,3) / ( k ♦ pow(sigma,2) ^ pow(epsilon,2) * 
pow(Tbase,9)); 

length = .82 * (-.0315113 * taper_ratio + 1.0425113) * Heat / ( sigma * epsilon * 
pow(Tbase,4)); 
db = length / 20016; 

length = length - db; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 


121 


if (depth < 10*base_width) Continue = No; 
if (length < 10*base_width) Continue = No; 

while (Continue == YesX 

length = length + db; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 

failures = 0; 
number_of_passes++; 

computet); 

while (qb == 0 && Continue == Yes){ 
failures++; 

length = length + db / 2; 

base_width = 2 Ap / ( length * ( taper_ratio + 1)) 
edge_width = taper_ratio * base_width; 
computet); 
if (failures == 3)C 
Continue = No; 
number_of_passes = 1;}} 


dq = qb - dq; 

cprintf("\r\nwidth = %4.3f length = %4.6f heat = %4.6f dq 
length, qb, dq); 
dq = qb; 


if (qb_max < qb )( 

max_length = length; 
qb_max = qb;) 

else 

Continue = No; ) 

if (number_of_passes > 2)C 
length = max_length; 

base_width = 2 * Ap / ( length * ( taper_ratio + 1)); 
edge_width = taper_ratio * base_width; 
computet); 

cprintf("\r\nWriting Rec # %d d = %4.3f b = %4.6f qb = %4 
j, base_width, length, qb, Tbase, Ap); 
out = fopen("OPTTRAP.MAT",‘*a‘'); 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper_rati( 
fclosetout);) 

elseC 

cprintf("\r\nIgnoring This Set Of Parameters"); 

return 0;} 

void compute(void){ 

double m, hr, n, ua, ub. A, B, C, D, fin_taper, I1UA, KOUB, lOUB, 

double qa, b1, dH, a, a_width, b_width, Ttip, temp, segment_diff, 

double b, inc, Tav, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 

int overflow, possible_overflow = No, number_of_tries = 0; 

b1 = length / ( 1 - taper_ratio ); 
fin_taper = atan2( base_width, 2 * b1 ); 
dH = length / N; 

K2 = 2 * epsilon * sigma * y * pow(Tbase,4); 

Tave = Tbase * pow(V, 0.25) + 0.0001; 

if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase • (Tbase * sqrt( 

else 


= %4.6f", base_width*1000. 


3f Tb = %4.1f Ap = %4.1f", 


%e, %e, %e\r\n", j, 

), base_width, edge_width ); 


K1UA, KOUA, lOUA, K1UB; 
HUB; 

100000; 


2 * V )))/2); 
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Ttip = lave + 0.0001; 
inc = .5 * Ttip; i=1; 

while ( fabsCoverall_diff) > ( Tbase * error) |] i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 

qa = sigma * edge_width * epsilon * pow(Ta, 4) * depth; 

a = b1 - length; 

while (i <= N && overflow == No){ 
segmented!ff = 10; 

a_width = edge_width * a / ( b1 * length); 
b^width = edge_width * ( a + dH ) / ( b1 - length); 

while (f absC segmented!ff ) > ( Tb * error) / N && overflow == No){ 

Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
n = hr / (k * sin(fin_taper)); 
ua = 2 * sqrtCa * n ); 
ub = 2 * sqrt( (a + dH) * n ); 

IlUA = Il(ua); 

KOUB = KO(ub); 
lOUB = lO(ub); 

K1UA = Kl(ua); 

KOUA = KO(ua); 
lOUA = lO(ua); 

KlUB = Kl(ub); 

HUB = HCub); 

A = ua * (!1UA*K0UB I0UB*KlUA); 

B = ub*ub*( lOUB’^KOUA - IOUA^KOUB) / ( 2*k*b_width*depth*n); 
C = 2’^k*a_width*depth*n*ub*(IlUB*KlUA - 11UA*K1UB)/(ua); 

D = ub*( 10UA*K1UB + 11UB*K0UA); 
temp =A*Ta+B*qa; 
qb = C * Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase + 50 |j Tb < 0)( 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;}} 

if (overflow == No){ 

Ta = Tb; qa = qb; a = a + dH; -*-+i;>} 


if (number_of_tries > 29)< 

T a=0;Tb=0;qb=0;return;) 

if (overflow == No){ 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overaIl_diff >0) Ttip=Ttip- inc; 
else Ttip = Ttip + inc;} 

else 

Ttip = Ttip - overall_diff; 
number^of_t ries++; 
if (Ttip < Tave) Ttip = Tave;}} 


return;} 
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j it ★★*★***★★★*★***********★**★***********★★****★★★★★★**★★★**★*★★****★★**★****★ 


PROFILE AREA Analysis Program Tb Determine The Optimum Geometry 
Longitudinal Fin of Triangular Profile 

ititititifkitititititititititititititicititititititititititititicitititiciticicititititititititititicititiciticiciciciticicicicic'kicic'kititititititic it j 


include <time.h> 
#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <conio.h> 
#include "bessel.h" 


#define Yes 1 

#define No 0 

^define N 99 

#define sigma 5.6696le'8 


double 

double 

float 

int 

void 


length, depth, edge_width, base_width, max_length; 

Ta, Tb, qb=0; 

V, qb_max, db, epsilon, k, error=0.01, taper_ratio, Tbase, <2, Heat, Ap; 
i. Exit, Continue, number_of_passes, j, failures; 
computeC); 




MAIN PROGRAM 




int rnainCX 
float dq; 

FILE *out; 
random! 2 e(); 
error = error / 100; 
cIrscr(); 

for (j=1; j < 102; j++)L 

number_of_passes = 0; 

epsilon = (float) ( random(98) ); 

epsilon = .01 + epsilon / 100; 

depth = (float) (20 + random(200)); 

depth = depth / 100; 

k = 10 + random(200); 

Tbase = 100 + random(500); 

Heat = (float) (25 + random(1000)); 

Heat = Heat / depth; 

V = (float) ( random(80)); 

V = .01 + V / 100; 

V = 0; 

taper_ratio = (float) (random(98)); 
taper_ratio = 0.01 + taper_ratio / 100; 
taper_ratio = 0; 

Exit = No; 

Continue = Yes; 
qb_max = -1; 
qb = 0; 
dq = 0; 

Ap = 1.304 * pow(Heat,3) / ( k * pow(sigma,2) * pow(epsiI on,2) * pow(Tbase,9)); 

length = .95 * Heat / ( sigma * epsilon * pow(Tbase,4)); 

db = length / 20016; 

length = length * db; 

base_width = 2 * Ap / length; 
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if (depth < 10*base_width) Continue = No; 
if (length < 10*base_width) Continue = No; 

while (Continue == Yes)( 

length = length + db; 
base_width = 2 * Ap / length; 

number_of_passes++; 
failures = 0; 

computet); 

while (qb == 0 && Continue == Yes)C 
failures++; 

length = length + db / 2; 
base_width = 2 * Ap / length; 
computet); 
i f (fai lures == 3)<; 

Continue = No; 
number_of_passes = 1;>} 


dq = qb - dq; 

cprintf("\r\ndelta = %4.3f length = %4.6f heat = %4.6f dq = %4.6f", base_width*1000, 
length, qb, dq); 
dq = qb; 


if (qb_max < qb ){ 

max_length = length; 
qb_max = qb;} 

else 

Continue = No; } 

if (number_of_passes > 2)C 
length = max_length; 
base_width = 2 * Ap / length; 
computet); 

cprintf ("\r\nWrit ing Rec %d d = %4.3f b = %4.6f qb = %4.3f Tb = %4.1f Ap = %4.1f", 
j, base_width, length, qb, Tbase, Ap); 
out = fopenC'OPTTRI.MAT","a"); 
edge_width = 0; 

fprintf(out, "%d, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e, %e\r\n*', j, 

depth, length, Ap, k, Tbase, epsilon, V, qb, taper_ratio, base_width, edge_width ); 
fclose(out);} 

elseC 

cprintf("\r\nlgnoring This Set Of Parameters"); 

j=j-i;» 

return 0;} 

void compute(void){ 


double 

double 

double 

double 

int 


m, hr, A, B, C, D, IIUA, KOUB, lOUB, KlUA, KOUA, lOUA, K1UB, HUB; 

fin_taper, b1, db, a, ua, ub, a_width, b_width; 

inc, Tave, overall_diff = 100, Tmax = 0, Tbv_min = 100000; 

Tav, Ttip, temp, qa, segment_diff; 

overflow, possible_overflow = No, number_of_tries = 0; 


edge_width = base_width / 100; 

db = length / ( N + 1 ); 

bl = length; 

length = length - db; 

fin_taper = atan2( base_width, 2*bl ); 


K2 = 2 * epsilon * sigma * v * pow(Tbase,A); 
Tave = Tbase * pow(V, 0.25) + 0.0001; 


if (V==0) 

Ttip = .8625 * Tbase * sqrt( 2 * V ) + ( ( Tbase - (Tbase * sqrt( 2 * V )))/2); 
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else 

Ttip = lave + 0.0001; 
inc = .5 * Ttip; i=l; 

while (fabsCoveral l_di f f) > ( Tbase error) | ] i <= N ){ 

overflow = No; i = 1; Ta = Ttip; Tb = Ta; 
segnient_dif f = 10; 

while (fabs(segment_diff) > ( Tb error) / N && overflow == No){ 

Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * powCTav, 3) - (K2 / (2 * Tav)); 
m = hr / (k * sin(fin_taper)); 
ub = 2 * sqrt( db * m ); 

lOUB = lO(ub); 

K1UB = KKub); 

HUB = IKub); 

C = 2'^k'*edge_width'* depth*pow(m,2)’*I 1UB/(ub * lOUB); 


temp = Ta / lOUB; 
qb = C * temp; 

segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 

if (Tb > Tbase +50 ]| Tb < 0)( 

overflow = Yes; possible_overflow = Yes; 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 

if (overflow == No){ 

Ta = Tb; qa = qb; ++i;} 

a = b1 - length; 

while (i <= N && overflow == No){ 
segment_diff = 10; 

a_width = edge_width a / ( b1 - length); 
b_width = edge_width "^ ( a + db ) / ( b1 - length); 

while (fabs(segment_diff) > ( Tb * error) / N && overflow == No){ 

Tav = (Tb + Ta) / 2; 

hr = sigma * epsilon * pow(Tav, 3) - (K2 / (2 * Tav)); 
m= hr/(k* sin(fin_taper)); 
ua = 2 * sqrt(a * m ); 
ub = 2 * sqrt( (a + db) m ); 

I1UA = Il(ua); 

KOUB = KO(ub); 
lOUB = lO(ub); 

K1UA = Kl(ua); 

KOUA = KO(ua); 
lOUA = lO(ua); 

K1UB = Kl(ub); 

HUB = I1(ub); 

A = ua * (IIUA’^KOUB + I0UB*K1UA); 

B = ub^ub^dOUB’^KOUA - I0UA*K0UB) / ( 2*k*b_width*depth*m) 
C = 2*k*a width’^depth^m^ub^dlUB’^KlUA - 11UA’^KlUB)/(ua); 

D = ub’‘( lOUA^KlUB + IIUB^KOUA); 

temp =A*Ta+B* qa; 
qb = C Ta + D * qa; 
segment_diff = ( temp - Tb) / 2; 

Tb = temp + segment_diff; 
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if (Tb > Tbase + 50 ]| Tb < 0){ 

overflow = Yes; possible_overflow = Yes 
if (Ttip < Tbv_min ) Tbv_min = Ttip; 
inc = (Tbv_min - Tmax) / 2; 

Ttip = fabs(Tbv_min - inc); 
if (Ttip < Tave) Ttip = Tave;>> 

if (overflow == No){ 

Ta = Tb; qa = qb; a = a + db; ++i;}> 

if (number_of_tries > 29)( 

Ta=0; Tb=0; qb=0; return;) 


if (overflow == No)< 

if (Tmax < Ttip) 

if (Tb < Tbase) Tmax = Ttip; 

overall_diff = (Tb - Tbase) / 2; 

if ( possible_overflow == Yes){ 
inc = inc / 2; 

if (overall_diff > 0 ) Ttip = Ttip * inc 
else Ttip = Ttip + inc;> 

else 

Ttip = Ttip " overalI_diff; 
number_of_tries++; 
if (Ttip < Tave) Ttip = Tave;)) 

length = length + db; 
return;) 
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APPENDIX C MATLAB EVALUATION M FILES 




MATLAB EVALUATION M FILE To Determine The Optimum Geometry 
Longitudinal Fin of Rectangular Profile 

load rect; 

D = rect; 

N = 100; 

sigma = 5.66961e‘8; 


L 

= D(; 

:,2); 

b 

= D(: 

,3); 

A 

= DC 


k 

= D(: 

:,5); 

To 

= D(: 

:,6); 

e 

= D(: 

,7); 

V 

= DC 

:,8); 

qo 

= D(: 

:,9); 

TR 

= D(: 

: JO) 

db 

= DC 

; J1) 

de 

= D(: 

: J2) 

I 

= [1: 

:101] 


C3 = A.* k.* (((sigma*e)."2).* To."9).* (L./qo).^3; 
Cl = b.* sigma.* e.* (To."4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (To.^5).* (L./qo).^2; 

dispCRECT V = O') 
disp(' ') 

disp('Constant One: Polynomial p is:') 

p = polyfit(TR, Cl, 1) 

fit1 = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum((Cl-fit1)."2)/(N-1)) 

pause 

dispC ') 

disp('Constant Two: Polynomial p is:') 

p = polyfit(TR, C2, 1) 

fit2 = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum((C2-fit2).^2)/(N-1)) 

pause 

dispC ') 

disp('Constant Three: Polynomial p is:') 

p = polyfit(TR, C3, 1) 

fit3 = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum((C3‘fit3).^2)/(N-1)) 
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I* ********lk**********************************W***** 

MATLAB EVALUATION M FILE To Determine The Opti 
Longitudinal Fin of Trapezoidal Profile 

•k-k-k-k-k-k-ifk-k-k-k-k-k-k-k'k-k'k'k'k'k'k'k'k-tt'k-k-k'k-tfk-k-k-tfk-k-k'k-k-k'k-tftfk-k-k-ifk-k-kifk 


load trap; 

D = trap; 

N = 100; 

sigma = 5.66961e-8; 


L 

= 

0(:,2); 

b 

= 

D(: 

,3); 

A 

= 

D(: 


k 

= 

0<:,5); 

To 

= 

DC 

,6); 

e 

= 

DC 

Jy. 

V 

= 

DC 

,8); 

qo 

= 

DC 

,9); 

TR 

= 

DC 

,10) 

db 

= 

D(: 

,11) 

de 

= 

DC 

,12) 

I 

= 

[1:101] 


C3 = A.* k.* (((sigma*e)."2).* To."9).* (L./qo)."3; 
Cl = b.’^ sigma.* e.* (To."4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (To."5).* (L./qo)."2; 

dispCTRAP V = O') 
dispC ') 

disp('Constant One: Polynominal p is:') 

p = polyfit(TR, Cl, 1) 

fiti = polyval(p,TR); 

dispCMean Square Error:') 

std^err = sqrt(sum((Cl-fit1).^2)/(N-1)) 

pause 

dispC ') 

disp('Constant Two: Polynominal p is:') 

p = polyfit(TR, C2, 2) 

fit2 = polyval(p,TR); 

dispCMean Square Error:') 

std_err = sqrt(sum((C2-fit2).^2)/(N-1)) 

pause 

disp(' ') 

disp('Constant Three: Polynominal p is:') 

p = polyfit(TR, C3, 1) 

fit3 = polyval(p,TR); 

dispCMean Square Error:') 

std_err = sqrt(sum((C3-fit3).^2)/(N-1)) 

pause 


mum Geometry 


************************ it j 
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j-k kk'kkkk'kk'kkk'k'k-kkkkkkkk'kkkkkkkkkkkkkk-k-kick-kk-kkkkkkkkk'k'kkkkkkkkkk-k-kkkk-kkkkk-kkk-k 


MATLAB EVALUATION M FILE To Determine The Opti 
Longitudinal Fin of Triangular Profile 


load tri; 

D = tri; 

N = 100; 

sigma = 5.66961e'8; 


L 

= DC 

:.2) 

b 

= DC 

:.3) 

A 

= DC 

:.«) 

k 

= DC 

.5) 

To 

= DC 

:,6) 

e 

= D(; 

:,7) 

V 

= DC 

:.8) 

qo 

= D(: 

:,9) 

TR 

= DC 

:,10) 

db 

= D(: 

:,11) 

de 

= D(; 

:,12) 

I 

= [1: 

101] 


C3 = A.* k.* (((sigma*e)."2).* To."9).* (L./qo)."3; 
Cl = b.* sigma.* e.* (To.^4).* (L./qo); 

C2 = db.* k.* sigma.* e.* (To."'5).* (L./qo). "2; 

dispCTRI V = O') 
disp(' ') 

disp('Constant One: Polynominal p is:') 

p = polyfit(TR, Cl, 1) 

fiti = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum((Cl-fi11).^2)/(N-1)) 

pause 

dispC *) 

disp('Constant Two: Polynominal p is:') 

p = polyfit(TR, C2, 1) 

fit2 = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum( (C2-f i t2).''2)/(N-1)) 

pause 

disp(' *) 

disp('Constant Three: Polynominal p is:') 

p = polyfit(TR, C3, 1) 

fit3 = polyval(p,TR); 

disp('Mean Square Error:') 

std_err = sqrt(sum((C3-fit3).^2)/(N-1)) 


mum Geometry 


kkkkkkifkifk-kifkkidckifkidikick kj 
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